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Abstract 



In this thesis we deal with the specific collective phenomena in condensed matter 
— striped-structures formation. Such structures are observed in different branches 
of condensed matter physics, like surface physics or physics of high-temperature su- 
perconductors. These quasi-one-dimensional objects appear in theoretical analyses 
as well as in computer simulations of different theoretical models. Here, the main 
topic of interest is the stability of striped structures in certain quantum models, 
where a tendency towards crystallization competes with a tendency towards phase 
separation, and some basic properties of these structures. 

We consider two strongly correlated two-component quantum systems, consist- 
ing of quantum mobile particles and immobile ones. By immobile particles we 
mean those particles whose occupation numbers are invariant with respect to the 
hamiltonian evolution. But this does not mean that these occupation numbers 
serve as a kind of a fixed disorder field. They are varied to reach the minimum 
of the energy of the total system (an analog of annealing). The both systems are 
described by Falicov-Kimball-like Hamiltonians on a square lattice, extended by 
direct short-range interactions between the immobile particles, that favor phase 
separation. In the first system the mobile particles are spinless fermions while in 
the second one they are hardcore bosons. We construct rigorously ground-state 
phase diagrams of the both systems in the strong-coupling regime and at half 
filling. Two main conclusions are drawn. Firstly, short-range interactions in quan- 
tum gases are sufficient for the appearance of charge stripe-ordered phases. When 
they occur, a first order phase transition between a checkerboard- crystal phase and 
phase-separated state (a segregated phase) is impossible: by varying the intensity 
of a direct nearest-neighbor interaction between the immobile particles, the both 
systems can be driven from a segregated phase via striped phases (for instance, 
via a diagonal-striped phase in the case of fermions, and via vertical- (horizontal-) 
striped phases in the case of hardcore bosons) to the checkerboard phase. Sec- 
ondly, the phase diagrams of the two systems (mobile fermions or mobile hardcore 
bosons) are definitely different. However, if the strongest effective interaction in 
the fermionic case gets frustrated gently, then the phase diagram becomes similar 
to that of the bosonic case. 

We show that any anisotropy of nearest-neighbor hopping eliminates the 7r/2- 
rotation degeneracy of the so called dimeric and axial-stripe phases and orients 
them in the direction of a weaker hopping. Moreover, due to the same anisotropy 
the obtained phase diagrams of fermions show a tendency to become similar to 
those of hardcore bosons. 

Finally, introducing a next-nearest-neighbor hopping, small enough not to de- 
stroy the striped structure, we examine rigorously how the presence of the next- 
nearest-neighbor hopping anisotropy reduces the 7r/2-rotation degeneracy of the 
diagonal-striped phase. The effect appears to be similar to that in the case of 
anisotropy of the nearest-neighbor hopping: the stripes are oriented in the direc- 
tion of the weaker next-nearest-neighbor hopping. 
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1 Introduction 



In the passing decade, some specific phases, having quasi-one-dimensional structure, the 
so called striped phases, have been a highly debated subject in condensed-matter physics. 
Apparently, a broad interest in such phases was initiated by reports presenting experi- 
mental evidence for the existence of charge stripes in doped layered perovskites, some of 
which constitute materials exhibiting high-temperature superconductivity P1I21- 

However, striped phases had been observed much earlier, for instance in physisorbed 
monolayers on metallic surfaces (see jS] and references therein), or in ultrathin magnetic 
films 13] . Many other instances of experimental observations of stripe-ordered phases are 
listed in |5^, where Monte-Carlo studies of formation of striped phases in a continuous 
gas with hardcore and short-range repulsive interactions are reported. Theoretical de- 
scriptions, including computer simulations, of these phenomena involve various kinds of 
classical lattice-gas models (also termed classical spin models) with competing interac- 
tions: for physisorbed monolayers see Ref. j^, for ultrathin magnetic films see Ref. |7j. 
Striped phases have been also studied in the framework of quantum-spin models, like 
XY model for instance fH]. Interesting theoretical considerations of stripe phases in a 
two-dimensional electron gas, physically realizable in MOSFET's, based on classical spin 
models with competing interactions, can be found in Ref. |9j. 

Quite interestingly, theoretical studies of striped phases in systems of strongly corre- 
lated electrons fOl El El Ej have preceded experimental observations of such phases. 
But only after those observations, the discussion became much more vigorous, and the 
nature of stripe-ordered phases started to attract attention of numerous researchers. In 
the context of the Hubbard model, a comprehensive review of the problem can be found 
in pi]. The existence of the same kind of striped phases was investigated also in the t-J 
model ^ni QEl EI- A bird's eye view on the problem of stripe-ordered phases in high- 
temperature superconductors, but emphasizing its general relevance for contemporary 
condensed- matter physics can be found in [IQ]. The general relevance of striped phases 
is underlined also in [I^, where they are viewed as an emergent phenomenon resulting 
from collective motions of microscopic particles, somewhat analogous to quasi-particles 
like phonons. 

The Hubbard or t- J-like models belong to the most realistic models, in the framework 
of which the problem of striped-ordered phases in doped layered perovskites can be 
investigated (see a review paper by Oles [El)- In the both models, the spin and the charge 
degrees of freedom are taken into account, and it is believed that it is the competition 
between these degrees of freedom that is decisive for the formation of striped phases. The 
major difference between Hubbard, t-J, and similar quantum-particle models for stripe 
formation, on one side, and the classical or quantum spin models (mentioned above) on 
the second side, is the role of the kinetic energy. It is absent in the second group of 
models while it plays a crucial role in the first group. The model studied by us in this 
thesis belong to the first group. 

The question of formation of striped phases is closely related to the question of relative 
stability of these phases against mixtures of an electron-reach phase and a hole-reach 
phase (segregated phases). Due to the tiny energy differences between both phases, the 
results obtained by means of approximate methods, which introduce hardly-controllable 
errors, are disparate. Therefore, as pointed out in [THl I2DI further careful studies are 
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necessary to settle the problem of formation of stripe-ordered phases. One of possible 
ways of attacking this problem is to formulate analog problems in less realistic but simpler 
models, where some control over the results obtained by means of various methods of 
statistical physics and many-body physics can be gained. This leads hopefully to a deeper 
insight into the, mentioned above, stability problem. 

Such an approach has been adopted by many researchers, using models of the both 
groups mentioned above. Here are two examples concerning classical spin models. To 
investigate the phenomenon of frustrated phase separation in high-temperature supercon- 
ductors. Low et al j2I] consider a spin 1 two-dimensional Ising model with short-range 
ferromagnetic coupling competing with long-range antiferromagnetic Coulomb interac- 
tions, where spin variables are coarsegrained representations of the local density of mo- 
bile holes. By means of combined analytical and numerical techniques they show that 
the transition between ferromagnetic and antiferromagnetic states proceeds via numer- 
ous phases, among which there are striped phases. In a recent paper, Valdez-Balderas 
and Stroud investigate, by means of Monte Carlo techniques, a competition between 
superconductivity and other types of order in two dimensions. Each site of the under- 
lying lattice, occupied by a classical XY spin (plane rotator), is interpreted by them as 
a mobile, positively charged superconducting domain, and the orientation of the spin as 
a phase of the superconducting order parameter of this domain. Vacant sites represent 
negatively charged nonsuperconducting domains. The nearest neighbor spins are cou- 
pled ferromagnetically (a representation of Josephson tunneling) and the occupied sites 
repel each other with a kind of a screened Coulomb interaction. By varying the relative 
strength of the two competing interactions, the system is driven from a phase separated 
state to a checkerboard-like state via a series of complex patterns of self-organization, 
including some kinds of striped phases. 

And here are examples concerning quantum-particle models. Buhler et al [22] have 
found, by means of Monte Carlo simulations, that upon hole doping antiferromagnetic 
spin domains and charge stripes, whose properties are in very good agreement with 
experiments, appear in a spin-fermion model for cuprates. Using the so called restricted 
phase diagrams, the stability problem of charge-stripe phases has been studied in the 
spinless Falicov-Kimball model by Lemahski et al [201 121] • In their study the formation of 
charge stripe-ordered phases can be looked upon as a way the system interpolates between 
a periodic charge- density wave phase (the chessboard phase) and the segregated phase (a 
mixture of completely filled and completely empty phases), as the degree of doping varies. 
A considerable reduction of the Hilbert space dimension in a spinless fermion model with 
infinite nearest-neighbor repulsion (as compared to a Hubbard model) has been exploited 
by Zhang and Henley [121 123 to study carefully, by means of an exact diagonalization 
technique, the formation of charge stripe-ordered phases upon doping. They addressed 
also an interesting question of the role of quantum statistics in the problem of striped 
phases, by replacing fermion particles with hardcore-boson ones. 

Our work has been inspired mainly by the recent studies of Lemahski et al [2011211 and 
by Zhang and Henley Pnil2n|- To investigate the key question, whether charge stripe- 
ordered phases are stable compared to a phase-separated state, we study rigorously the 
strong-coupling limit of an extended spinless Falicov-Kimball model on a square lattice, 
with mobile particles being spinless fermions or hardcore bosons. The usual spinless 
Falicov-Kimball Hamiltonian (such as that studied in [20] ) has been augmented by a 
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direct, Ising-like interaction between the immobile particles. 

In the framework of striped-structures formations, one of the interesting questions, 
the influence of hopping anisotropy on striped phases, was investigated by means of the 
Hartree-Fock method in |2Z]. They showed that under an anisotropy in hopping 
intensities, stripes become oriented in the direction of weaker hopping. This problem is 
investigated here rigorously, using the introduced model. 

The thesis is organized as follows. In the next section we give some insight into 
the Falicov-Kimball model, providing its main properties, and we introduce technics 
which are used for deriving our results. Then, in Section 3 we present the considered 
models and construct phase diagrams due to truncated effective interactions. After that, 
in Sections 4 and 5, we analyze the influence of hopping anisotropy on striped phases. 
Finally, we draw conclusions and provide a summary. Various technical details a placed 
in Appendices (from A to E). Our original results are contained in Sections 3,4, and 5, 
and have been published in [2H1 1211 IHO] , respectively. 
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2 The Falicov— Kimball models 



In this section we give the definition of the spinless Fahcov-Kimball model and it basic 
properties. First of all, we discuss possible interpretations of the model in respect to 
history of its development. We also shortly mention about today's topics of interest. 
In the second paragraph we present some useful properties of the model. The following 
paragraphs deal with technics we use to investigate the model: we make a sketch of a 
perturbation scheme and discuss the m-potential method in more details. 

2.1 Definition 

The spinless Falicov-Kimball model, on arbitrary lattice A is described by the Hamilto- 
nian: 



where sites on lattice A are denoted as x,y, . . ., whose number is |A|, and the first sum- 
mation is over the pairs of different sites x and y, with each pair counted once. The 
system described by Hamiltonian ((T)) consists of two sorts of spinless particles: quantum 
which can hop (usually between nearest-neighbors sites) on the lattice with hopping am- 
plitude txy = t*^ (we call them electrons) and classical which are immobile (ions). The 
electrons, as quantum particles, are described by creation and annihilation operators 
on site x: and Ca,, respectively, which satisfy the appropriate commutation relations: 
canonical anticommutation relations for spinless fermions or commutation relations of 
spin 1/2 operators 5*+, S~ , for hard-core bosons. The total electron- number operator 
is A^e = Tlx^x = Tlix^x^x and corresponding electron density is pe = Ne/\A\. Classical 
ions, are described by pseudo-spin (or simply, spin) Sx on cite x: if Sx = 1, site x is 
occupied by an ion and, if Sx = —1 the site x is empty. Ions are described in terms of 
spins for convenience: instead of spin variable Sx we can use occupation-number variable 
^x — (■5a + l)/2. The collection of spins values {sx}x£A is called the ion configuration. The 
total number of ions is Ni = ^^.n^ = '^xi^^ + aiid the ion density is pi = Ni/\A\. 
There is no direct interaction either between the electrons or between the ions, although 
there is the simple on-site interaction between the two kinds of particles with coupling 
constant U. It is clear that the particle-number operators A'e, Ni, and spins Sx are 
conserved. Therefore, the description of the classical subsystem in terms of the ion con- 
figurations 5* = {sx}xeA remains valid. Whenever periodic configurations of pseudo-spins 
are considered, it is assumed that A is sufficiently large, so that it accommodates an 
integer number of elementary cells. 

The Falicov-Kimball model was first considered by Hubbard [31, and Gutzwiller (S2] 
as a simplification of the one-band Hubbard model. 



where txy = 0, and txy = txy for arbitrary sign of coupling constant U . This is an 
approximation of the Hubbard model in the case of infinitely heavy spin-up electrons. 






(2) 



o-=T,J, {x,y)ck 
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In 1969 this model was introduced by Falicov and Kimball in order to explain metal- 
insulator transitions in mixed valence compounds of rare-earth materials jnS|- Experi- 
ments suggested that these transitions occur only due to the interactions between elec- 
trons: localized immobile /-electrons (ions) and d-electrons in Bloch-like states. The 
model described by dH) is the simplification of the original model, suggested in 

In 1986, Kennedy and Lieb reinvestigated Falicov-Kimball model as a "static Hubbard 
model". It was used in [Slj and jSHI to study crystallization. Some rigorous results 
concerning the model were presented in these works for the first time. Mainly, they 
proved the existence of the phase transition in the model for any value of U and for 
sufficiently low temperature. 

In the last years the Falicov-Kimball model is intensively studied in the frame of the 
dynamical mean field theory, where properties of the system could be calculated exactly 
in the limit of infinite dimension (see review [^). 

We would like to mention about the numerical studies of the Falicov-Kimball-like 
models, which are provided on finite clusters and are performed by Farkasovsky and 
co-workers (see for e.g. jSZ], [SH], EHl)- 

Due to the relative simplicity of the model, a variety of modifications have been 
applied to it. First of all it can be considered on different types of lattices. In jlUj, the 
model was considered on a triangular lattice. Due to increasing interest of such systems 
on lattices such as Kagome or pyrochlore, studies of the Falicov-Kimball model on them 
are going to be instructive. 

An interesting problem arises on changing the statistic of quantum particles. In jlH! 
the electron subsystem, whose particles obey Fermi statistic where replaced by a system 
of hard-core bosons. Since two bosons cannot occupy the same site, because of hard-core 
interaction, the subsystem consists of interacting particles. A comparison of this two 
statistics might be useful for computer simulations, where fermions are often replaced by 
hard-core bosons, in order to make computations easier. 

One of the alterations, which was studied recently, is the so-called correlated hopping, 
i.e. the hopping intensity t^y depends on occupation numbers of ions on sites x and y. 
Influence of the correlated hopping on the phase diagram of the simple Falicov-Kimball 
model was investigated in ^I] in the case of ID model, and in jl21 and jlSI in 2D. 

Another simple modiflcation, which perhaps makes the model more realistic, is in- 
troduction of a spin degree of freedom. It was studied by Brandt et al. in [HI and 
jj45j. Recently, using perturbation technique, the Falicov-Kimball model in case of infl- 
nite Coulomb repulsion between classical ions with different spin orientation at the same 
site was investigated in jlHl- The model with an additional spin-depended Ising-type 
interaction (Ising-Falicov-Kimball model) was studied recently in |47j . 

Historical review and list of new trends, as well as extended comments, concerning 
the fleld of the Falicov-Kimball model are given in a short guide jlH] . 

2.2 Symmetries and other basic properties 

We are interested in ground-state properties of the Falicov-Kimball model. We describe 
its ground state from the point of ion subsystem, so the main question is the ground- 
state conflguration of ions. Here the model on a square lattice is considered, though many 
properties of the model can be extended to a wider class of latices. The most convenient 
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ensemble for our purposes is the grand-canonical ensemble. Let, 

H {lie, f^i) = HpK - fJ'eNe - fJ^iNi, (3) 

where chemical potentials of electrons and ions, respectively. 

First, we notice, that for a fixed configuration {sx}x(^a Hamiltonian (^, in the case 
when electrons obey Fermi statistic, describes a system of free Fermi particles in an 
external field U ■ Sx- The Hamiltonian (^J is the second-quantized version of the single 
particle Hamiltonian 

h = -T + US, (4) 

where h is |A| x |A|-matrix with the elements hxy = txy + UsxSxy However, the configura- 
tion of ions {sx}x£A is not fixed in the system: the configuration of ions should minimize 
the whole Hamiltonian H [n^, Hi), and the main question of describing the ion subsystem 
is to find such a configuration for specific values of (/ie;/^j)- This follows from the form 
of a partition function, 

Za = ^ Tr exp [-PH (/ie, /ij)], (5) 
s 

where the sum runs over all the configurations and the trace is taken over the fermion 
Fock space (/3 is the inverse temperature). In other words, if Es {fj,e, fJ^i) is the ground- 
state energy of H (/ig, /ij) for a given configuration S of the ions, then we are looking for 
the ground-state energy of H (/ig, /ij), which is given by. 

Eg (/ie, /ii) = inin Es (/ig, /li) . (6) 

The minimum is attained at the set G of the ground-state configurations of ions. The 
task is to determine the set G in the (/ig, /ii)-plane. 

For this purpose, it is useful to examine the spectrum of the Hamiltonian (in the 
case of fermions). When max \txy\ = t and \U\ > zt, where z is the coordination number 
of the underlaying lattice {z = 4, for square lattice), there is a gap in the spectrum 
{—\U\ + zt, \ U\ — zt) (the so-called universal gap). This gap is of great importance, 
because if we suppose the chemical potential of electrons fj,e to be in the universal gap, 

fieeHU\ + ^t,\U\~^t), i\U\>zt), (7) 

then the ground-state energy expansion, we use below, is absolutely convergent, and 
uniformly in A. Quite interestingly, a similar condition holds in the case when quantum 
particles are hard-core bosons, 

^Me^ {-\U\+4zt,\U\-4zt), {\U\>4zt). (8) 

This condition was obtained by Messager and Miracle-Sole in 02] , using a cluster expan- 
sion. Under conditions ((Tj) and (jHI), for U positive (negative) and in the zero-temperature 
limit, the half-filling condition, Ng + Ni = |A| {neutrality condition, Ne = Ni), is implied, 
either for fermions or hard-core bosons. 
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In studies of grand-canonical phase diagrams an important role is played by uni- 
tary transformations {hole-particle transformations) that exchange particles and holes: 
c^Cx — 1 — c+Ca;, for electrons, and — >■ —s^, for ions. 

The peculiarity of the model is that the case of attraction {U < 0) and the case of 
repulsion {U > 0) are related by a unitary transformation (the hole-particle transfor- 
mation for ions): if 5* = {sx}a;gA ^ ground-state configuration at {fie, fJ'i) for U > 0, 
then —5* = {— 5^;}^^^ is the ground-state configuration at {fie, ~f^i) for f/ < 0. This 
property holds for both statistics, and does not depend on the form of hopping term. 
Consequently, without any loss of generality one can fix the sign of the coupling constant 
U . In the following paragraphs we choose f/ > and express all the other parameters of 
the Hamiltonian Q in the units of f/, i.e. formally we set f/ =1, preserving previous 
notations. 

In the case, when only nearest-neighbor (n.n.) hopping terms are present, i.e. t^y = 
for all {x, y)i, ^ ^ 2, this transformation for some (/ig, leave the Hamiltonian H {fie, f^i) 
invariant. An example of such transformation for the electrons is given by —>■ exCx, 
where ex = 1 for bosons; for fermions e^; = 1 at the even sublattice of A and e^; = — 1 at the 
odd one. Since Hpx is invariant under the joint hole-particle transformation of mobile 
and localized particles (when only n.n. hopping is non-zero), H {fi^, fii) is hole-particle 
invariant at the point (0,0). At the hole-particle symmetry point, the system under 
consideration has very special properties, which simplify studies of its phase diagram 
|34j . Moreover, by means of the defined above hole-particle transformations one can 
determine a number of symmetries of the grand-canonical phase diagram jHO!- With the 
sign of U fixed, applying joint hole-particle transformation for electrons and ions, one 
concludes that there is an inversion symmetry of the grand-canonical phase diagram. If 
S' is a ground-state configuration at {fie,fii), then —S is the ground-state configuration 
at {—fte,—fii)- Therefore, it is enough to determine the phase diagram in a half-plane 
specified by fixing the sign of one of the chemical potentials. 

Let us note, that in the case when the next-nearest-neighbor (n.n.n.) hopping ampli- 
tude is nonzero, txy ^ for {x,y)2, Hamiltonian Q is not invariant with respect to the 
joint hole-particle transformation for any values of fie and /ij. 

Let us consider symmetry properties of the model with the following restrictions: 
n.n. hopping is supposed to be independent of direction, i.e. txy = t, for {x,y)i, and 
n.n.n. hopping amplitude differs txy = in two different directions. Moreover, we 

consider the system of fermions only. Then applying the hole-particle transformation for 
electrons, X, and the hole-particle transformation for ions, one finds that 

if 5* is the ground-state configuration at {t,t^,t^, fie, fii), then —S is the ground-state 
configuration at {t, —t^, —t^, —fte, —fii) ■ Thus, it is enough to consider only one sign 
of n.n.n. hopping. We shall consider the case of positive n.n.n. hopping intensities: 
> 0. Using another hole-particle transformation for electrons, Cx — > c^, and the 
hole-particle transformation for ions, we obtain that if S is the ground-state configuration 
at {t, t+, t-, fie, fii), then —S is the ground-state configuration at {—t, — 1+, — t-, —fie, —l^i)- 
Applying consecutively the two joint (with respect to electrons and ions) hole-particle 
transformations, we obtain that if 5* is the ground-state configuration at (t, t+, t_, fie, fJ^i), 
it is also the ground-state configuration at {—t,t+,t_, fie, fii). So the relative sign of n.n. 
and n.n.n. hopping amplitudes does not play any role in the model. 

Above, we have briefly discussed those properties of our model, which shall be useful 
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in the sequel. For more extended picture of the Fahcov-Kimball model, its modifications, 
results and applications, readers are referred to a series of reviews jlHI, [HI], [S2] and jHSl- 

2.3 Strong-coupling expansion of the ground-state energy at 
half-filling — the method 

There are different ways of obtaining ground-state energy expansion. In the case of 
fermions, when zt < 1 and |/ie| < 1 — zt, the ground-state energy of the model can be 
written in the form [201 1^ '■ 

Es {fJ^e,fJ^i) = -^Tr|/l| - (/ii - fJ.e)Ni - /ie|A|, (9) 

where the one-particle Hamiltonian h = —T + US is given by (0)), and \h\ = ^/h?. The 
trace of \h\, can be rewritten in the form 

Ti\h\ = TrV/ + A, (10) 

where, I is the unity operator, A = T'^ — {TS+ST) stands for small perturbation, S'^ = I. 
Expanding the square root into a power series of A, one obtains the ground-state energy 
expansion. Such a "square root expansion" was applied in jHO]- 

An alternative way for fermions, the "resolvent expansion", was used in jlHI- Ground- 
state energy Q can be rewritten in the form 

-{fli- fie)N,- fi,\A\, (11) 

where C is a contour in the complex plain enclosing all the negative eigenvalues of h. In 
order to calculate the above integral trace of resolvent is expanded into a power series of 

txy 

Unfortunately, the above methods are not applicable in the case when quantum hop- 
ping particles are hard-core bosons. To obtain the ground-state energy expansion in this 
case, a closed-loop expansion is used, which was first applied in ^\ for fermions, and 
then used for hard-core bosons in |40j . 

We use another, more powerful method, the method of unitary-equivalent interactions, 
which is suitable not only for both systems, fermions and hard-core bosons, but can also 
be applied to a wide class of Hamiltonians. It was proposed in a series of papers [54j, 
|55j . [Sn] and summarized in [2Z|. The main advantage of this method is that the unitary 
transformations are applied locally, and that is why the method is well-defined: the 
expansion is convergent, uniformly in |A|. 

The idea of the method is the following. Considering the Hamiltonians which are sums 
of local terms, H = Hx {X is a finite set of sites in A), they are supposed to consist 
of two parts: classical, which is diagonal in some basis, and small quantum perturbation, 
which is off-diagonal. With A being the small parameter considered in the model, we are 
looking for some unitary operator Uf^\X), which makes the original Hamiltonian i?A(A) 
block-diagonal to order A" in some basis, i.e. 

iff (A) = f/f (A)/fA(A) f/f * (A) = <(A) + vt\X), 



Es (/Lte,/Ui) 



Sx + 



dz 
27ri 



Tr 
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where the leading block- diagonal part Hq^{X) contains terms up to the n-th order in A, 
and off-diagonal part V^"'^(A) consists of {n + l)-th and higher order terms. Following 
this procedure, the convergence of the expansion can be rigorously proven. Moreover, by 
applying quantum Pirogov-Sinai theory, it was established that some parts of the grand- 
canonical zero-temperature phase diagram of the leading part is modified "slightly" when 
the whole expansion is taken into account at zero or low-temperature regimes [SZj. 

Here we present a scheme of obtaining the ground-state energy expansion for the 
Falicov-Kimball model up to the fourth order. We consider the case where nearest- 
neighbor and next-nearest-neighbor hopping intensities are anisotropic. 

2.3.1 Local projections 

First we consider the system, where quantum hopping particles obey Fermi statistic. 
Following inZ], we rewrite FK Hamiltonian Q in the form convenient for perturbation 
(up to some constant), i.e. "classical part" plus "small quantum perturbation": 

Hfk = Ho + Vu (12) 

where the classical part, 

Ho = J2^0x, (13) 

with the classical on-site potential, 

$0x = 2<<. - /ie<. - funi; (14) 
and the quantum perturbation, 

{x,y)i,hCA (a;,j/)i,„cA {x,y)2,+Ck {x,y)2,^CA 

i {x,y)^cA 

Operators Q(x,y)i are defined in the following way: 

Q{x,y)i,h ~ {^x^y + CyCx) , 

Q(x,y)i,v ~ i^x ^y ~^ ^y ^x) ) 

Q {x,y)2,+ ^+ ipx ^y ~^ ^y ^x) ; 

Q{x,y)2,- = t- [c+Cy + C+Cx] , (16) 

and the role of the small parameter A is played by t = max {t/i, t+, (the so-called 
strong- coupling regime) . The method is applied in the situation where for certain value of 
energy of the classical part there are infinitely many states. Then, it is possible that the 
quantum perturbation reduces this degeneracy by sphtting the energy levels. Let us look 
at the classical part, or the "zeroth order" . The phase diagram of the classical interaction 
(Fig. is partitioned into four parts. On each boundary line between any two phases 
the degeneracy is infinite. We are interested in the case of half-filling which corresponds 
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{1e, 1i} 


{Oe, 1i} 




(2,2) 




{1e, Oi} 


(0.0) 






{Oe, Oi} 







Figure 1: Ground-state phase diagram for Hq = Yl, {'^'>^t:'>^x ~ f^en^ — /^i?^^)- Numbers of 

electrons and ions at each site of lattice A, which determine the state of the system, are 
given in curly brackets. 



to the boundary line /i,. = (0 < /ie < 2), the configurations on this line with the 
lowest energy are characterized by the condition: each site of lattice A is occupied by one 
particle only, either ion or electron. 

On the line /le — we introduce the projection onto the ground state at site x: 

Px = < (1 - <) + < (1 - <) = « - <)' • (17) 

The projection onto excited states ~ 1 — P^. It is easy to extend these projections to 
the ground-state projection on y C A: 

In this case unity is given by ly = Py -|- Py, and hence the definition of the projection 
onto exited states at the set Y: 

P^^ty- po (note, that i^y 7^ H (1^) 

xeY 

For calculation we need the following simple relations: 

c+po = pic+ 

^x-^ x X^X 1 

ctPl = P°c+, (20) 

U y ^ X, then c+ or commute with any projection at site y. These relations give 
immediately very useful properties for any operator Q{x,y)i, with pair x y on lattice A: 

P{x,y)iQ{x,y)i = Q{x,y)iPxPy^ 

pLhQi^^vh = Q(.,.>, (1 - ^X) ■ (21) 
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Another pairs of relations which correspond to ()20|) and ()21|) are obtained by Hermitian 
conjugation. 

Using these properties and the partition of unity, for Y = {x,y)i, we can rewrite 
Q{x,y)i- 

Q(x,y), = Q^x,y), + Qtl,y), + <5(x,y>,5 (22) 

for i = {Ih, If, 2+, 2—}; where 

Q^x,y), = P(x,y),Q{x,y),P(x,y), = 

Qtx,y). = P(x,y),Q{x,y),P(x,y), + P{x,y),Q {x,y)rP{x,y)i 
a; J/ x-'- y ' x-'- y {x ,y} i J- x-'- y t 

Q(is/>, = P{x,y),Q{x,y),Plx,y),- (23) 

So the whole Hamiltonian (|12j) is rewritten in the form: 

Hfk = Ho + + V'\ (24) 

where, 

{x,y)i.h {x,y)i,v {x,y)2,+ {x,y)2,- 

« {x,y)i 

{x,y)i,h {x,y)i,v {x,y)2,+ (x,y)2,- 

= Y.Y.Q"U- (25) 

2.3.2 First-order unitary transformation for fermions 

We are looking for a unitary transformation, that eliminates the first-order off-diagonal 
term V^^. The transformation can be written in the form: 

f/(i) = exp(S«), (26) 

where the operator 5*^^^ = S'^^'> {th,ty,t+,t-) is of order one, i.e. it depends on t^t^t^^tl 
with a + b + c + d= 1. Applying the first unitary transformation and using the 
Lie-Schwinger series, the transformed Hamiltonian can be written as 

= Hfk +[S'^'\ Hfk] + [S^'\ [S'^'\ Hfk] ]+.•• = 

oo ^ 

= E;^^d"5«(^^^)' (27) 

n=0 
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where, the operation adA{B) is defined as follows 

ad°A(5) = 5, 
ad^v4(5) = [A,B], 
ad"A(5) = [A,a(r-'A{B)]. 

Rewriting fTf\ . 

i/d) = Ho + 1^°^ + + ad^^^i) (ifo) + ad^S^i) + V') + 

+ y' ad"5(^) (i/o + + V') , (28) 



n=2 



and choosing the operator S^^^ in the way to eliminate V^^, 
we obtain: 



n! \n + 1 / 



n=2 

In the above formula, 



Now our problem is to find the operator S^^\ If we suppose it to be a sum of local 
operators, 

in order to satisfy the condition [Hq, S''-^-'] = V^^, ^t^^y)- '^^^ chosen in the form: 

^{x,yh = ad"^i^o,(x,j/>, {Q^{x,y),) , (30) 



where ifo,(x,2/)i = '^'oa; + Introducing the spectral decomposition of ifo,(a;,i/>i, 

j 

S^^^ can be rewritten as 

p(i) r)01 p(fc) 

.(1) \^ \^ ^{x,y)y{x,y),-^{x,y). 



q(l) _ ST^ ST^ {x,y)i^ {x,y), {x,y)i 
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In the spectral decomposition of $01 there are three eigensubspaces with eigenvalues 
Ex^ and with the corresponding projections Px \ / = 0, 1, 2: 



Ei'^ = 0, 

= 2(l-/x), 



P(0) 

X 

P(l) 

X 



< + <-2«, 



where fi = fii = fie- Note, that Pi^"* = P°, and Pi^'' + Pi'''' = P^. Then, the spectral 
decomposition of Ho^(^x,y)i consists of six eigenvalues e'^^^^^^ and eigenprojections Pf^y)-^ 
1 = 5: 



>(2) 



I P(l) P(0) 

^ X J/ 5 

+ Pi'^P^'\ 



^(0) 

{x,y)i 


= -2/x. 


p(0) 


- P(0) p(0 

x ?/ 


^(1) 

{x,y)i 


= -/^, 


p(l) 


- P(0)p(l 
X ^ y 




= 2-3^, 


p(2) 


- P(0)p(2 


E?^ , 
{x,y)i 


= 0, 


p(3) 

{x,y)i 


_ P(l) P(l 
^ X ^ y 


{x,y)i 


= 2(l-/i), 


p(4) 


- P(l) P(2 
X y 


Ef , 

{x,y)i 


= 4(l-/i). 


p(5) 


- pmp{-i 

X y 



I P(2) p(l) 
' X y 1 



Considering the products P\i]y)P%,y)^- ^x,y),, 



Pf'^ \ , one finds that only two of them are nonzero: 



P 



(0) 



{x,y)P%.,y)E^^,y)^ ^ud P'^ly)P%,^y)P'^ly)^. for any i. Using the projections defined above. 



>(o) 



we can write the expression for operator S**^^-*: 

S^^^ = ~9 5Z X] {ExPyQ%,y),PlPl - PlPlQ%,y),PxPy) 



(32) 



To proceed, we have to consider commutators ad^y4(P) = AB — BA, where A = 
Yj{x,y), '^{x,y)^ and B = Yjj Y.{m,n)j B{m,n)j- First, let us note that the commutator for 
\_Ai^x,y)ii B(Jn,n)J^^ vauishcs unless {x,y)i fl {m,n)j 7^ 0. In this case, the product. 



AB 



x:(x,y)i 



{x,y)iB[x,y)i 



x:{x,y)i z:{y,z)j 



The first term in the brackets runs over all the possible pairs of sites, i.e. i = {Ih, Iv, 2+, 2—}. 
In the second term we introduce the extended indices of summation, x : {x,y)i and 
z : {y, z)j, in order to stress that first sum runs over all a; G A and pair (x, y) is the n.n. 
pair of the i-type, and the second sum runs over z (shown on Figs. I2l3li|) . such that pair 



{y, z) is the n.n. pair of j-th type. Whenever we skip the index of direction (like '/i' or 
'— ') this means that the summation is over all the possible directions. So, the notation 
J2{xy)2 stands for sum over all n.n.n. pairs, each counted once. 

The second-order term is V2, which after somewhat lengthy calculations can be written 
in the form: 



9 Ex Py Q{x,y)iPx Pi 



i x:{x,y)i 
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Figure 2: The set of paths over which the sum ^^..^^.^^^ J2z-{yz)i ^^^^ (ESI)- Sums run 
over the set of paths shown here, and their translations. Sites x and z are at the ends of 
each path (x can be chosen arbitrary), and site y is always between them. 
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Figure 3: The set of paths over which the sum J2x-{xy)i 'l2z-{yz)2 ^^^^ (jHHj) . Sums run 
over the set of paths shown here, and their translations. Sites x and z are at the ends of 
each path {x can be chosen arbitrary), and site y is always between them. 



i^j x:{x,y), z:{y,z)j 

^2"=E(^ E PxPlQlx,y),PlPij- 
i^j \ x:{x,y)i 

~2 ^ ^ ] |^-Px-Py-P2;'5(a;,j/)iQ(y,2)j -Pr-Py-Pz ~ Px Py Pz Q (x,y)iQ (y,z)j Px Py Pz ~ 

x:{x,y)i z:{y,z)j 

~ Px Py Pz Q {y,z) j Q {x,y)i Px Py Pz ~^ Px Py Pz Q {y,z) j Q {x,y)i Px Py Pz ~^ 



PxPyPz [Qix,yhQiy,z),+Qiy,z),Qix,y).] PXP'- 



- PxPyPz [Q{x,y),Q{y,z), + Q {y,z),Q {x,y),] PxP^Pz 



We have got all the elements which are necessary for obtaining the effective Hamilto- 
nian up to the third order (the second-order ground-state energy expansion): 
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Figure 4: The set of paths over which the sum ^^..^^.^^^ '^z-{yz)2 ^^^^ Sums run 

over the set of paths shown here, and their translations. Sites x and z are at the ends of 
each path {x can be chosen arbitrary), and site y is always between them. 



Hence, 



i x:{x,y)i 



In the ground state + = 1, so 



QUy). = t\ [nl (1 - nl) + (1 - <)] = - 1) , 



{a;,2/>i 



and finally the second-order ground-state energy expansion is of the form: 

4 



(2) 



1 



{x,y)i 



where /i = /ij — /Xe and = /Xj + /ig. 



2.3.3 Second-order unitary transformation for fermions 

Applying the second unitary transformation, we obtain: 



(34) 



(35) 



(36) 
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where V^^^ is the remainder of the fifth order: 



\n + l J 

+ ad^5(2) i. ad"5(^) (;^^°' + ^'') ) + 

2! ^ 2 2 2^^, 1^^^^ 

In order to ehminate the second-order ofi^-diagonal term V^^ ^ we put: 



yOl ^ yll 



then 
and 



ii-(2) = Ho + V^^ + 1/2°° + 1^2" + V3 + Vi + V^^\ 



(37) 



The operator 5"^^^ can be chosen in the form: 

c(2) _ V V V 9^^) 

x:{x,y)i z:{y,z)j 

where, analogously as for 5"*^^^ 



01 



{x,y)iU{y,z)j 2-^ {x,y)iU{y,z)j -^{a) — E^^^ (^.2/>iU(2/,z>j ' 

<^,P {x,y)iU(y,z)j (x,y)iU(y,z)j 



The spectral decomposition used in the above formula, is given by 

^ox + ^ov + ^oz = y^ Ej""^ \ ,,/ \ -P/°\ ,,/ s ■ 

ua; I U2/ I vz £^ {x,y)i\J{y,z)j (x,y}iU{y,z}j 
a 

Operator 5"*^^^ is constructed in the similar way as was. It reads: 

x:{x,y)i z:{y,z)j 

-P(x,y)^U(y,z)j [Q{x,y)iQ{y,z)j + Q {y,z)jQ {x,y)i] P{x,y)My,z)^ ■ (39) 
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The next step, consists in calculating double commutators, [A, [B, C]] that appear 
in in (j37j) and (|38|) . Each commutant X, of these commutators, has the form X = 
'^x-{x y)i -^{x,y)i- Applyiiig the same argument as before (commutators vanish unless com- 
mutants are defined on subsets intersection of which is not empty), we restrict ourselves 
to the products where the sets of sites, on which each multiplier is defined, forms a con- 
nected path, i.e. each site on a path has nearest neighbor, or next-nearest neighbor, 
or both. Another thing, worth mentioning, is that, even considering the fourth order 
only, we do not need to calculate all the expressions in (jHTj) and (jHHj) . The result, we 
are interested in, is projected onto ground state, hence, some of the terms vanish. The 
following property of operators Q(x,y), is useful. Because of relations it is clear, that 
if we consider the path, path = {ii,i2, ■ ■ ■ , iq}, the expression 

PpathQ(il,i2)aQ(i3,ii)b ' ' ' Q {iq-i,iq) z^path 7^ 0) (40) 

only if each of indices ij in the subscript of Q repeated even number of times. With this 
condition satisfied, let us note that any operator = P^AP^ + P'^AP^ or = PMP^ 
is equal to zero after projecting onto ground states. We will use these properties during 
calculations third- and fourth order terms. 



Third order. The calculation of the third order terms is the calculation of V3. First 
note, that in our case there is only one type of path that consists of three pairs of n.n. or 
n.n.n sites and product of three operators Q on it satisfies the condition (pOj) (see Fig. E)). 
In addition, after ground-state projection only V^^" is nonzero. Using the multiplication 



ii O i ! > {(- 



Figure 5: The set of paths of Pi-type which appear in the third-order calculations, 
rule in a way described above, we obtain: 

^ ~ T ^ ^ Px Py iQlxyQlyzQ2xz ~l" QlyzQlxyQ2xz ~\~ QlxyQ2xzQlyz~\~ 



Pi 



QlyzQ2xzQlxy ~l~ Q2xzQlxyQlyz ~l~ Q2xzQlyzQlxy) ^x^y^z ' 



or in terms of spin variables: 



yOO 
^3 



Pi,- 



(41) 



where the set of paths Pi splits into two parts, Pi^+ and Pi,-, where n.n.n. part of the 
path has the slope +1 and —1, respectively. The sp^ . stands for the product of spins 
which belong to Pi,j. This is the third-order correction to the effective Hamiltonian. 
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Fourth order. To consider fourth-order products, we start with describing connected 
paths that satisfy condition (^01) for a product of four operators Q. There are three types 
of such paths: X-type path, i.e the path consisting of the pair (n.n. or n.n.n.) of two sites 
X; XF-type path, i.e. the path consisting of three sites that formed by two pairs X and 
Y; XYZW-tjpe path, i.e. the path consisting of four sites which are in the form of closed 
loop and are formed by four different pairs of sites, — X,Y, Z and W. The corresponding 
product of Q's to the X-type path is Q\. The product proportional to Q^Qy corresponds 
to XF-type path. The product of the type QxQyQzQw corresponds to the XYZW- 
type path. Only these types of products contribute to the ground-state projection of V4. 
Another thing connected with calculations, is that some terms, like ad^S'^^'*(V2'^^), vanish 
on projecting onto the ground-state. Performing all necessary calculations, we finally 
obtain: 

a x:{x,y)a 

+ X] X] {^Q'lx,y)cQ\y,z)p- Q{x,y)c.Q'ly,z)i3Q{x,y)c.- Q{y,z)pQ\x,y)^Q{y,z)^ - 

a,/3 x:{x,y)a z:{y,z)p 

^ [l^ - [Q{x,y)^Q{z,j)^ + Q{z,j)^Q{x,y)^] " [Q {y ,z) {x ,j) ^ + Q {x ,j) {y ,z) + 
i=2 x-.xePi \ 

+ 2 ■ [Q{y,z)^Q(x,j)^ + Q(x,j)^Q(y,z)^] ■ [Q(x,y)o.Q(z,j)^ + Q (z,j)^Q (x,y)c,] + 

+ [Q{x,y)M{y,z)^ + Q{y,z)f>Q{x,y)^] " [Q {z,j)^Q {x,j)^ + Q {x,j)^Q {z,j)^] + 

+ [Q{y,z)pQ{z,j)-, + Q(z,j)-,Q(y,z)p] " [Q {x,j)^Q (x,y)^ + Q (x,y)^Q (x,j)^] + 

+ [Q{z,j)^Q{x,j)^ + Q{x,j)^Q{z,j)^] ■ [Q{x,y)^Q{y,z)p + Q {y,z) {x,y)o.] + 

+ [Q(x,j)^Q(x,y)c, + Q(x,y)o.Q(x,j)^] " [Q (y,z) f^Q (z,j)-, + Q (z,j)^Q (y,z) ■ 

The paths Pi are shown in Fig. IHl We group the terms and rewrite it using spins variables. 
Note, that during this procedure the summation over Pg is replaced. Summing up, the 
fourth-order correction to the ground-state energy expansion, reads 
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Figure 6: The sets of paths P2, ■■■ , Pq which appear in the fourth order calculations. 



+ I —thl - — 



' 16 ^ 16 



{^■,y)3,v 



{x,y)4.h- 



{^■,y)4,v-] 



{'■c,y)4,v- 



+ 3^^+^- E (5^^3 + 1) + ^^^^+^- E + 1) + 3^^'^+^- E (5^^*- + 1) + 
+ yg^'^I E (5^^5,.+ + 1) + E (5^^5,.- + 1) + 



The whole ground-state energy expansion up to the fourth order, 



(42) 



4^^) = Ho + ^2°° + V^^ + + , (43) 

where V^^, V^^ and V^*^" are given by ()34|) . ()41|) and (021), respectively. The remainder 
i?*^^) is independent of the chemical potentials, and collects all the terms proportional to 
tlt\,t'+ti , with a + 6 + c + = 5, 6, . . .. 



2.3.4 The case of hard-core bosons 

On replacing the hopping fermions with hard-core bosons, the commutation rules change: 
operators of creation and annihilation at different sites commute. Using the same technic 
as in the case of fermions we obtain the ground-state energy expansion in the case of hard- 
core bosons. For the system of hard-core boson we consider only the case t+ = t_ = 0. 
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The second order correction to the "classical" part Hq remains the same as for fermions. 
The third order disappears (due to the absence of n.n.n. hopping, on the square lattice 
the path that satisfy condition (j4m) for a product of three Q operators does not exist), 
and the fourth-order correction, reads: 



2.3.5 Building the ground-state phase diagram 

We make some remarks on the expression of the ground-state energy expansion (j43j) 
in both cases: hopping fermions and hard-core bosons (in the case of hard-core bosons 
V^^ = and V^^ in is replaced with {V^^)hcb, given by (HH)). In the above expressions 
we keep even the terms, that are independent of spin variables, hence have no influence 
on the phase diagram and so can be omitted. The effective Hamiltonian ()43|) depends 
only on /i = /ij — /ie, which we call the unique chemical potential parameter. 

The procedure of building the phase diagram is recursive. The phase diagram of the 
effective interaction truncated at the order k' is constructed on the basis of the phase 
diagram obtained at the preceding order k: the conditions imposed on the ground-state 
configurations by the k-th order effective Hamiltonian have to be obeyed by the ground- 
state configurations of the A;'-th order effective Hamiltonian. In other words, the k'-th 
order terms of the expansion cannot change the hierarchy of configuration's energies 
established by the k-th order effective Hamiltonian; they can only split the energies of 
configurations in cases of degeneracy. 

2.4 m-potential method 

To proceed, we need to introduce some definitions. First, as was already mentioned, 
we denote the configuration of spins on the lattice A as = {sx}xeA, and takes two 
values, either +1 or —1, at any site x. We call Sa the restriction of the configuration 
S to the subset of sites A C A with 1^41 sites, i.e. Sa = {sx}xeA- Instead of talking 
about functions of spin variables s, it is convenient to use the function of configuration 
S (spin Sx depends on the configuration at site x, Sx). For A C A, let /a be a function 
of configurations such that. 



2.4.1 Definition of m-potentials 

The concept of the m-potential was introduced in [SH], and its idea is the following. 
Suppose that for some A C A, it is possible to write a Hamiltonian, which is a function 
of configurations, in the form. 




5^ ^xSy + 



(44) 



fA{S) = fAiSA), 



for any S. 




(45) 



A 
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where by the sum we understand the summation over all translations of A in A. 

Assume, that there exists a certain set of periodic configurations {S} (global configura- 
tions) which minimize the potential Ha{S) for {Sa\ (local configurations). If it holds for 
all translations of A in A, then we say that Ha is an m-potential and the set {S} is the 
set of the ground-state configurations. 

Let us present two examples. Consider the antiferromagnetic Ising model on a square 
lattice (J > 0), 

H I sing J ^ ^ SxSy. (46) 

The potential H(^x,y)i{S) = JsxSy is an m-potential. Indeed, the local minimum for this 
potential is for configurations of nearest-neighbor pairs of spins where one is up and the 
second is down: { — 1,1} or {1,-1}. There exist two global configurations, called the 
chessboard configurations, which minimize the potentials for all {x,y)i C A: if 

we choose any nearest-neighbor pair of sites, then the local configuration of the pair 
will be one of those which minimize the H(^x,y)i- In other words, if we can extend the 
local configuration (which is the ground-state configuration locally) to the global one, 
in a way that other local configurations with higher energy do not appear, we call the 
local potential Ha the m-potential, and the global configuration is the ground-state 
configuration because every potential under the sum has minimum on it. 

Now, we consider the same Ising model pUj) on a triangular lattice. As in the case of 
square lattice, local nearest-neighbor configurations remain the same. But the problem 
appears while trying to extend the local configurations to the whole lattice A: there is 
no such global configuration, which is ground-state configuration locally. This means 
that -f^(x,y)i is not an m-potential in this case (rewriting (jmj) as a sum over elementary 
triangles, it appears that the corresponding potential is m-potential, however the number 
of ground states is infinite). 

Practically, potentials depend on some parameters and often they are m-potentials 
only for some range of these parameters. In order to find the complete phase diagram, 
i.e. all ground-state configurations for the whole range of the parameters, we are forced 
to build new m-potential. There is no recipe for constructing m-potentials. The only way 
to deal with the problem, as far as we know, is to introduce the so-called zero-potentials. 

2.4.2 Definition of zero-potentials 

The idea is simple: if some local potential does not satisfy the m-potential condition, then 
let's change it by adding zero, which means that we change it locally, but in a way that 
the total Hamiltonian remains the same. So the most general definition of zero-potential 
reads: a potential Ka-, is a zero potential if it satisfies the zero-potential condition: 

Ka{S) = 0, for any S. (47) 

A 

Now, it is clear that adding such a potential to the Hj>^, 
HA{S) = J2HA{S) = J2iHA{S) + KA{S)) = J2HA{S), for any 5, (48) 

A A A 

modifies only the local potential Ha- Under the zero-potential condition (j47|) . Ka has 
no definite structure. In the sequel we shall consider two ways of constructing Ka- 
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2.4.3 First kind zero-potentials 



In physics literature one can find two different ways of constructing zero-potentials. The 
first method was proposed in [SUl, where a part of the ground-state phase diagram of the 
Falicov-Kimball model was obtained for the first time in the fourth order. To illustrate 
the idea of this method we provide an example. Let us consider Hamiltonian with a 
local potential Ha which is not an m-potential. Now we take some function (or several 
functions) of configurations on a set of sites A. For example, one of the simplest choices 
is a ^ Sxi where a is a number. We add and subtract this function, so the Hamiltonian 

does not change. 



if A = X j if A + a X Sa; - a X 



(49) 



Next we chose a new set of sites A' (or several sets), and rewrite the Hamiltonian in the 
form: 



ifA 



X ffA + a X s J - a ■ I A| X s^. = X ffA + a 5Z ) + 

A \ xeA J xgA A \ xeA / 

4/ \ I I ^ ^.c 4' A A' 



(50) 



The problem is to find such a value of the coefficient a that converts potentials Ha and 
if A' into m-potentials. Here it means that the m-potential condition must be satisfied 
for the both types of sets A and A': the restriction of a global ground-state configuration 
to any A and A' should minimize the corresponding potentials H^ and H^'- 



2.4.4 Second kind zero-potentials 



Another method, was proposed in [59] and in |40] . In contrast to the first method, where 
zero-potentials are of the form /a(5') — /a(5') = (i.e. the zero-potential condition is 
satisfied locally), no such restrictions are applied to its form. However, it is convenient to 
choose a zero-potential in a way that leaves the Hamiltonian's symmetry unchanged. We 
show here how to obtain zero potentials for the case when A is a 3 x 3-plaquette, later 
on called the T-plaquette (Fig. [7j). We assume that ifr does not change under rotations 



Figure 7: The sites of T-plaquette are labeled from the left to the right, starting at the 
bottom left corner and ending in the upper right one. 



by 7r/2 and reflections. 
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Any function f{S) can be rewritten in a form, 



f{S) = J2KS')xiS',S), (51) 

S' 

where S' runs over all possible configurations, A(5") is a value of function f{S) for specific 
configuration 5" and x{S', S) is an indicator function of configuration S', i.e. it is equal 
to 1 when S = S' and is zero otherwise. 

The construction of an indicator function x{S', S) for specific configuration 5" is 
simple. At a site x the indicator function for spin up is Xx{+'^, S^) = (1 + Sx)/2, for spin 
down is 1, S^) = (1 — s^)/2. The indicator function at A, is defined as 

Xa{S'a,Sa) = IIxx{S'x,Sx). (52) 

Performing the summation in (jM^ and taking into account the form of the indicator 
functions (j^ . any function /a(5'a), can be written as, 

Sz + ... + Xa'Wsx = 

xGA {x,y}<zA {x,y,z}<zA x&A 

= 5^AxSx(5), (53) 

XQA 

where Ax are the values which determine the function Ja- Notation sx{S) stands for 
the product of all spins which belong to X, sx = IIxgx -^"^^ X = 0, the coefficient 
Ax = Ao and sx = 1- 

Another way to show the identity (j53|) and give a formula for finding the coefficients 
in is to consider the functions sx{S) instead of indicator functions x{.S',S). If 
is the number of spins in Sx oriented down, then sx{S) = (— 1)^~. Let us define the 
scalar product in the set of functions of configurations: 

/i\ |A| 

{fA,9B)=l-) $^/a(5')^7b(5'), (54) 
where S' runs over all possible configurations on A. Now we show, that 

Let us consider the product, 

sxiS)sYiS) = sx.y{S), 

which follows simply from the form of sx{S); here X - Y = (X U F) \ (X fl Y). Using the 
fact, that for any set X, 

^ / 2l^l, ifX = 

l^sx{S) = ^ (56) 
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we obtain expression ()55|). Since the functions sx are orthonormal and there are as many 
sx as xa, they constitute an orthonormal basis, and any function can be expanded 
as follows, 

fAiS) = ^ XxSxiS) 

XCA 

which is equivalent to (|^. where 

\x = {sx,fA). (57) 

We are looking for a zero-potential Kt{S) in a form (jSHjl . We can always set Aq = 
although it can change the Hamiltonian only by a constant. So for a T-plaquette, which 
we will use later, we choose the zero-potential in a form: 

Kt{S) = Yl ^asa{S) = 

0^ACT 

= ^xSx + ^xySxSy + ^xyzSxSySz + . . . + XtSt- (58) 

3;6T {x,y}cT {x,y,z}cT 

Now applying to these objects the symmetry operations of the Hamiltonian and the 
zero-potential condition we can eliminate some of the arbitrary coefficients A^. For 
instance, we consider the first term: 

Ksx = Aisi + A2S2 + A3S3 + A4S4 + A5S5 + XqSg + A7S7 + A8S8 + A9S9. (59) 

x£T 

Demanding ^^^^ X^Sx to be invariant with respect to symmetry operations (rotations 
or refiections), we obtain 

Ai = A3 = A7 = Ag, A2 = A4 = Ae = Xs- (60) 

The zero-potential condition (jTTj) now reads: 

T xGT i x£A 

and holds for any configuration S. That directly implies the condition 

$^A. = 0. (61) 

i 

Expressing A5 in terms of Ai and A2, 

A5 = -4 (Ai + A2) , (62) 
only two coefficients remain independent, and the expression is of the form: 

A^s^. = Ai (si + S3 + S7 + sg - 455) + A2 (S2 + S4 + S6 + S8 - 4^5) . (63) 

xeT 

26 



Similarly, we can obtain other zero-potentials applying the procedure sketched above to 
products of nearest-neighbors spins, next-nearest-neighbors, etc. Finally, we can rewrite 
Kt{S) in the form 

KT{S) = J2a,4\S), (64) 

i 

where 

(S) = Si + S3 + S7 + Sg - 4S5, 

k!^\S) = S2 + S4 + se + ss - 4s5, (65) 

are the only two symmetric zero-potentials of the form (|59|) . For a T-plaquette, there are 
16 independent symmetrical zero potentials. The complete list of them is presented in 
Appendix A. 

Now, having zero potentials, we are trying to find such a vector a = {ai,a2, ■ ■ ■) 
that makes the total potential (the original potential Ha plus the zero potential Ka) 
an m-potential. Components are called the zero-potential coefficients. Both, in the 
first and the second methods, there are no restrictions for a but we suppose, that they 
are functions of the Hamiltonian's parameters Pi. For convenience, we write the set of 
parameters as a vector p = {pi,p2, ■ ■ ■). 

The both methods have their own advantages and disadvantages. For our purposes 
the second one is more suitable for two reasons. First of all, adding a zero-potential 
conserves the symmetry of the system. As a result, the number of block configurations 
whose energies we have to compare is reduced significantly. The second reason is that 
extending some local configurations to the whole lattice is much simpler in the case when 
we have one type of blocks, say T-plaquettes only. If other blocks are present, like A' in 
the first method, then we have to check if potential Ha' in (jSUjl . is an m-potential for 
a configuration we have built using ground-state configurations of the potential Ha, or 
vice versa. 

2.4.5 Some remarks on constructing zero-potentials 

Let the Hamiltonian, a function of configurations, be the linear function of parameters pi 
(by linear function we understand function which is convex and concave simultaneously), 
and it is possible to rewrite it in the following form 

Ha{p;S) = J2Ha{p-^S), (66) 

A 

where Ha is linear in p and is not an m-potential. Adding the zero-potential Ka in the 
form (jMjl . we suppose that a zero-potential vector a is the linear function of a vector of 
parameters p, 

H^ip; S) = J2 iHA{p\ S) + KA{a{p)\ S)) = Ha{p\ d{p)\ S), (67) 

A A 

and we are looking for a vector a for which the above potential is an m-potential for any 
vector of parameters p. Unfortunately, we hardly ever succeed in finding a zero-potential 
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vector which converts iJ^ into an m-potential for the whole range of the parameters; 
usually coefficients d{p) can be found in a certain range of the parameters only. So 
sometimes we need several sets of zero-potential coefficients to cover the whole range of 
the considered parameters. Such an approach is used in Section 3. 

In studies of the minimum of the potentials, the following observation is useful. Con- 
sider two points in the parameters space, x and y, and such a', a" in each point, respec- 
tively, that the potentials Ha{x] a'; S) and HA^y; a"; S) are minimized at some sets of 
local configurations Ga{x, a') and GA{y, a"), respectively, and Ga{x, a') nGA{y, a") 7^ 0- 
There exist only one linear function d{p) on the line segment between points x and y 
whose values are equal to a' at x and a" at y: 

d(tx + (1 - t)f) = t ■ a' + (1 - t) ■ a", for any t G [0, 1]. (68) 

This definition ensures that potential Ha is linear in p: 

Ha (tx + (1 - t)y, a (tx + (1 - t)y) ; = 
= t ■ Ha{x, 5) + (1 - t) ■ HA{y, d{y); S), for any t E [0, 1]. (69) 

We claim that a set of block configurations that realizes the minimum of energy 
Ha{p', d{p); S) on the line segment between x and y, is 

GA{tx + {l-t)y,d{tx + {l-t)y)) = GAix,d{x)) nGA{y,d{y)), 

for < t < 1. 

Indeed, let S E G'^(x, a(x)) fl GA{y,d{y)), i.e.: 

Ha{x, d{x); S) = min #a(^, d{x)] S) and ^^(y, d{y)] S) = min HA{y, d{y); S). (70) 
Then by and (fTO)). 

Ha (tx + (1 - t)y, d {tx + (1 - t)y) ; = 

= t ■ Ha{x, dix); S) + {l-t)- HA{y, d{y)- S) 
= t ■ min Ha{x, d{x); S) + {1 — t) ■ mm HA^y, d{y)] S) 

^ t ■ Ha{x, dix); S) + {l-t)- HA{y, d{y)- S) 
= Ha (tx + (1 - t)y, d {tx + (1 - t)y) ; , 

for any configuration S and any t E [0, 1], that is, 

Ga{x, d{x)) n GA{y, d{y)) C Ga [tx + (1 - t)f, d {tx + (1 - t)y)), 
for any t E [0, 1]. 

Now we shall demonstrate that for t E (0, 1) the opposite inclusion holds. Let 

S eGa [tx + (1 - t)y, d {tx + (1 - t)y)) 
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that is, 



Ha (tx + (1 - t)y, a {tx + (1 - t)y) ; ^ 



^ HA(tx + {l-t)y,a{tx+ {1 - t)y) ; S 



for any S. In particular, the above inequahty holds for S G GA{x,d{x)) fl GA{y,ci{y)), 
and can be rewritten as 



t 



Ha{x, d{x); S) — Ha{x, d{x) \ S) 



+ 



+ {l-t) HA{y,a{y)-S)-HA{y,d{y)-S) ^0, (71) 



for any t G (0, 1). Simultaneously, since each of the energy differences (in square brackets) 
in (ffTjl is nonnegative. 



Ha{x, d{x); S) — Ha{x, d{x) \ S) 



+ 



+ {l-t) HA{y,a{y)-S)-HA{y,a{y)-S) ^0. (72) 



By (HD) and (|72D, 



Ha{x, d{x); S) — Ha{x, d{x); S) 



+ 



+ (l-t) HA{y,diyy,S)-HAiy,d{y);S) 



0, 



with t > and 1 — t > 0, hence 

HA{x,d{x);S)-HA{x,d{x);S) = and HA{y,d{y}-,S) - HA{y,d{y)-,S) = 0, 
which implies in turn that S G Ga{x, d{x)) fl GA{y, d{y)). Therefore, 

GA{tx + {1 -t)y,a{tx + {1 -t)y)) C Ga(x, a(f)) n ^^(y, tt(y)). 

This proof has a simple geometrical interpretation (Fig. |H)). All energies are linear 
functions of t, so if there exist a point at t G (0, 1) where some block configuration S S 
has a lower energy than that of S, then at least at one of the points, x or y, the energy 
of S is lower than the energy of local ground state configuration. 

So, there is no need to find zero-potential coefficients everywhere: in order to find 
them for some interval (x, y) , it is enough to find them at x and at y, where the sets of 
local ground state configurations are G^ and Gjj, respectively. The construction described 
above guarantees that with a {tx + (1 — t)y) for t G (0, 1), the set of local ground state 
configurations Gp.y] on the line segment (x, y) is given by G[g^jfi = Gg fl Gg. Such an 
approach is used in Section 4 and 5, where zero-potential coefficients are given only on 
some semi-infinite lines and at certain points. 
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H,(s) 




1 t 

Figure 8: The illustration of the proof concerning a. 
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3 Stability of striped phases 



In this Section^ we consider two strongly correlated quantum systems, described by 
Falicov-Kimball-like Hamiltonians on a square lattice, extended by direct short-range 
interactions between the immobile particles. In the first system hopping particles are 
spinless fermions while in the second one they are hardcore bosons. Ground-state phase 
diagrams of the both systems in the strong-coupling regime and at half-filling are con- 
structed rigorously. Two main conclusions are drawn. Firstly, short-range interactions 
in quantum gases are sufficient for the appearance of charge stripe-ordered phases. By 
varying the intensity of a direct nearest-neighbor interaction between the immobile parti- 
cles, the both systems can be driven from a phase-separated state (the segregated phase) 
to a crystalline state (the chessboard phase) and these transitions occur necessarily via 
charge-stripe phases: via a diagonal striped phase in the case of fermions and via vertical 
(horizontal) striped phases in the case of hardcore bosons. Secondly, the phase diagrams 
of the two systems (mobile fermions or mobile hardcore bosons) are definitely different. 
However, if the strongest effective interaction in the fermionic case gets frustrated gently, 
then the phase diagram becomes similar to that of the bosonic case. 

3.1 Introducing the modifications to the Falicov— Kimball model 

The task is to generate the physical situation, where the striped phases appear, in the 
framework of Falicov-Kimball model. For technical reasons, our analysis is restricted to 
the case of half-filling, where the well-defined expansion can be written. However, it is 
now well established (proven) [26, £0, that, in the fermion spinless Falicov-Kimball model 
the phase-separated state, the segregated phase, is stable only off the half-filling. Thus, 
an additional interaction is needed to stabilize the segregated phase at half-filling. 

In the regime of singly occupied sites and for strong-coupling, the strongest effective 
interaction in the Hubbard model is the Heisenberg antiferromagnetic interaction. Doping 
of holes can be seen as means to weaken the strong tendency toward antiferromagnetic 
ordering and to make possible a phase separated state with hole-reach and electron-reach 
regions. Under similar conditions, in the spinless Falicov-Kimball model an analogous 
role is played by the Ising-like n.n. repulsive interaction between the immobile particles 
that favors a chessboard-like ordering. The effect of weakening of the strong tendency 
toward chessboard ordering can be achieved by an extra Ising-like n.n. interaction that 
compensates the strongest repulsive interaction, and consequently permits the system to 
reach a phase-separated state, the segregated phase. On varying, in a suitable interval 
of values, the corresponding interaction constant, which is our control parameter, we 
can study how the both systems "evolve" from the crystalline chessboard phase to the 
segregated one. 

Here, we would like to make a remark concerning physical interpretations of our ex- 
tended Falicov-Kimball model. The classical spin models, like those studied in flR 
can be considered as classical effective models of crystallization that contain the necessary 
"ingredients" of such models. Namely, in these models there are two competing inter- 
actions, one favoring phase separation into macroscopic regions, which are either empty 
or completely filled with particles, the other favoring periodic arrangements of particles. 

^This section is based on pS], 
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Our extended Falicov-Kimball model, with hopping fermions, can be thought of as a 
microscopic quantum model of crystallization. In this case, the necessary "ingredients" 
are: the tendency towards crystalline orders, coming from the kinetic energy of itinerant 
electrons and from the screened (in our case on site) electron-ion Coulomb interaction, 
a kind of Van der Walls attractive forces, favoring phase separated state, and finally the 
hard-core repulsion provided by the underlying lattice and the Pauli principle. 
In what follows, the modified Hamiltonian reads: 



H = HpK + V, 



(73) 



where Hpx is given by (P), and a direct Ising-like interaction V between the immobile 
particles is of the form 



V 



W 



E 



SxSy 



16 



E 



SxSy. 



(74) 



Here W stands for an intensity of direct n.n. interaction between ions. We also in- 
clude a subsidiary direct interaction between the immobile particles, an Ising-like n.n.n. 
interaction with intensity e, much weaker than the n.n. interaction. This interaction 
can reinforce or frustrate the n.n. interaction, depending on the sign of its interaction 
constant. It appears that on varying W, the systems (fermion or boson) characterized 
by different values of n.n.n. interactions, may undergo different sequences of transitions 
between the chessboard phase and the segregated one. We find, in particular, that if we 
set appropriate values of the n.n.n. interaction in the fermion system and in the boson 
one, then the both systems settle in the same phases, for typical values of the control 
parameter. On the other hand, the n.n.n. interaction enables us to demonstrate that the 
ground-state phase diagrams in the cases of mobile fermions and mobile hardcore bosons 
are definitely different. 

To simplify the model, we consider the case where only n.n. hopping is nonzero and 
does not depend on direction. So we put th = = t, and t+ = t_ = 0. We are interesting 
in investigation of the ground-state phase diagrams of our systems, so the grand-canonical 
formalism is used. We use the expansion (jlHj) which was obtained in the strong-coupling 
regime and at half-filling. In this case, with the expansion terms up to order four shown 
explicitly, it reads: 



fermion 



e _ 9t^ w 



3t^ _ e 

Te ~ 16 



+^ E '^'y + E (^'p + 1) + ^^'^ 



(75) 



P2 



for fermions, and 



E 



boson 
S 



E^- + 



4E 



4 16 



+ 



5t^ 
16 



e 

16 



{^,y)2 



SxSy 



16 



E(5 + 3p)+/2(^), 



(76) 
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for hardcore bosons, up to a term independent of the ion configuration and the chemical 
potentials. The remainders R^'^\ R^^\ are independent of the chemical potential /i and 
the parameters W and e, and collect those terms of the expansion that are proportional 
to t^"", with m = 3,4, .... 

Due to the convergence of the expansions (f73j) and (fTHj) . it is possible to establish 
rigorously a part of the phase diagram (that is the ground-state configurations of ions are 
determined everywhere in the (/ig, /ii)-plane, except some small regions), by determining 
the phase diagram of the expansion truncated at the order k, that is according to the fc-th 
order effective Hamiltonians (^^^'"™°"){'^') an^l ^j^hoson\^{k) order to construct a 

phase diagram according to a k-th order effective Hamiltonian, we use the m-potential 
method. 



3.2 Zeroth- and second-order effective interactions — phase di- 
agram 

As was discussed above, up to the second order the effective Hamiltonians for fermions 
and hardcore bosons are the same. Hence, the discussion that follows applies to the both 
cases, and the common effective Hamiltonians are denoted as Ef^ (/i), E^g^ (/i). 
In the zeroth order, 



^ {x,y)i {x,y)2 



where 



P2 



= J] (s, + 1) + ^ ^ s,sy -^Yl ^-^y^ 



'P2 



16 ^ " 16 

{x:,y)i (x,y)2 



and the primed sums in (|7Hjl are restricted to a plaquette P2- For any e > 0, the 
plaquette potentials Hp^ are minimized by the restrictions to a plaquette P2 of a few 
periodic configurations of ions on A. Namely, — the empty configuration (5*+ — 
the completely filled configuration), where Sj, = — 1 at every site {s^ = +1 at every 
site), and the chessboard configurations S^f^, where = e^, and 5'°^, (where = —^x), 
with e^. = 1 if X belongs to the even sublattice of A and = —1 otherwise. Moreover, 
out of the restrictions of those configurations to a plaquette P2, S^ip^, >S'+|P2, S^^p^, and 
S°^p^, only four ground-state configurations on A can be built, which coincide with the 
four configurations named above. Clearly, this is due to the direct n.n.n. attractive 
interaction between the ions. The section of the phase diagram in the (W, /i, e) space by 
a plane with e > 0, according to the effective Hamiltonian E^^^ (/x), is shown in Fig. El 
In the second order, 



e 

16 ^ -^"^ 

{x,y)i {x,y)2 



(79) 
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Figure 9: Ground-state phase diagram of 
E^g^ (/i) for fermion and hard-core boson sys- 
tems. 



Figure 10: Ground-state phase diagram of 
Eg (^) for fermion and hard-core boson sys- 
tems. 



where HpJ, 

and the primed sums are restricted to a plaquette P2- Since E^g^ differs from E^g^ by the 
strength of n.n. interactions, the only effect on the phase diagram is that the coexistence 
lines are translated by the vector (— 2t^,0,0). The section of the phase diagram in the 
(W, jjL, e) space by a plane with e > 0, according to the effective Hamiltonian e'^^ {jj), is 
shown in Fig. ^1 



H 



(2) 
P2 



s, + 1) + 



+ 



W 



8 16 



3.3 Fourth-order effective interactions — phase diagram 

It follows from the investigation of the phase diagram up to the second order that the 
degeneracy is finite, independently of the size of A, everywhere except the coexistence 
lines in the [e = 0)-plane, where it grows exponentially with |A|. Only there the effect of 
the fourth order interactions can be most significant. The meeting point of the coexistence 
lines in the [e = 0)-plane, where the coexistence line of 5*+ and sticks to the stability 
domain of chessboard configurations, appears to be particularly interesting. In that point, 
the energies of all the configurations are the same. In what follows, we shall study the 
phase diagrams of spinless fermions and spinless hardcore bosons up to the fourth order, 
in a neighborhood of radius O(t^) of the point (— 2i(:^, 0, 0). In this neighborhood, it is 
convenient to introduce new coordinates, {W,^,e) — > {uj,6,e), 

W = -2t^ + t^u, fi = t^6, e = t^e, (81) 

and a new (equivalent to t-independent effective Hamiltonian, 

^^ferrmonp (0,5) = ^ (^/™n^ (4) _ ^g2) 
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Then, in the spirit of the m-potential method, we express (^f/"^^'''"*°")(4) ^]^g poten- 
tials (i7/-™°-){4)^ 



where, in terms of the new variables, 

{Hl~p = -6 (.5 + 1) + ^ - ^) E" '^'y + ^ (3 - .) ^^sy + 

In flHH|) and (jHH), T stands for a T-plaquette of a square lattice (see Fig. [7j). The double- 
primed sums are restricted to a T-plaquette. For bosons, we introduce {H^°''"'Y^^ and 
in the same manner as for fermions, with 

(H'-T^ = -5 (.5 + 1) + ^ - E" '^'y + ^ (^ - E" ^^^^ + 

^ E" ^-^^^ - ^ E" (^^2 + 5) • (85) 

{x,y)3 P2 

We have to search for the lowest-energy configurations of (^^|'^^™°")(4) ^^^^^ ^j^bosonY-i) 
among all the configurations. Consequently, the potentials (^ff^^''™°")(4) g^j^d (Hlj^'^"^)^^^ 
have to be minimized over all the T-plaquette configurations. There are (up to the 
symmetries of Hq) 102 different T-plaquette configurations, shown in Fig.|2niin Appendix 
B. 

Unfortunately, in contrast to the lower-order cases, the potentials (^ff/'^^™°"y4) 
(^j^bosonyi) f^-yj^Yn out to bc the m-potentials only in a small part of (cj, 5, £:)-space. This 
difficulty can be overcome by introducing the zero-potentials k!^\ that are invariant 
with respect to the symmetries of H and are given by (jMj) . Coefficients at, depending 
on (a;,5, e) in general, have to be determined in the process of constructing the phase 
diagram, and the potentials are listed in Appendix A. We set all = for z ^ 6. 

Following (|iS|). new candidates for m-potentials can be introduced, 

(i7/™")W = Y (^(/7^™")(^) + K^^^^ , (86) 

T 

with an analogous representation of {H^^^°'^)^^\ and now the task is to minimize the po- 
tentials (i7^''")W + ir^''^ and (i7|?^°") ^ + JsT^^^ over all the T-plaquette configurations. 

In our study of the ground-state phase diagrams we limit ourselves to the {6 = 0)- 
plane, where the both considered systems are hole-particle invariant, and to the {e = 
0)-plane. Our analysis of the minima of the T-plaquette potentials, (^ff^'^''™°"y4) 

(if|P^°")(4) augmented by the zero-potentials shows that these planes are partitioned 
into a finite number of open domains each with its unique set of periodic ground- 
state configurations on A, denoted also by S^- There is a finite number, independent of 
|A|, of configurations in and they are related by the symmetries of H. The last two 
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statements do not apply to only one of the domains, denoted §^2, which will be described 
in the sequel. 

The domains §d are characterized as follows: at each point p {p = {uj,e) or p = {uj,6)) 
of a domain there exist a set of coefficients {ai{p)} such that the corresponding 
potentials are minimized by a set Std(p) of T-plaquette configurations. Moreover, from 
the configurations in Std{p) one can construct only the configurations in The set of 
the restrictions to T-plaquettes of the configurations from §£,, Sd\t, is contained in each 
set Std(p) with p E So- In Table ITD of Appendix E we mark by asterisk the cases, where 
the set Std{p) contains, besides S_d|t, some additional T-plaquette configurations. 

Specifically, for 5 = the ground-state phase diagrams due to the effective Hamilto- 
nians (^ffj^^™'*""'^^) ^^^^ j^^6oson^(4) ^^^^ shown in Fig. ITTl and Fig. IT^ respectively, while 

the corresponding ground-state phase diagrams for e = 0, in Fig. El and Fig. EI In 
Fig.Hniwe display the representatives of the sets §d of ground-state configurations. That 
is, the remaining configurations of Sd can be obtained easily by applying the symme- 
tries of Hq to the displayed configurations. The domain Sseg, having no representative 
in Fig. El consists of the two translation- invariant configurations S+ and S- , related by 
the hole-particle transformation. 

Only in the domain 8^2, which appears in the phase diagrams shown in Fig. El 
that is off the hole-particle symmetry plain, the situation is not that simple. In §^2 one 
can distinguish two classes, §d2a and 8^26, of periodic configurations with parallelogram 
elementary cells. A configuration in Sd2a consists of vertical (horizontal) dimers of filled 
sites that form a square lattice, where the sides of the elementary squares have the length 
2a/2 and the slope ibl. In a configuration of §^26, the elementary parallelograms formed by 
dimers have the sides of the length 2^/2 and the slope ±1, and the sides of the length ^/TO 
and the slope ±1/3. Two configurations, one from §d2a and one from §d2b, having the same 
kind of dimers (vertical or horizontal), can be merged together along a "defect line" of the 
slope ±1 (dashed line in Fig. EI), as shown in Fig. EI without increasing the energy. By 
introducing more defect lines one can construct many ground-state configurations whose 
number scales with the size of the lattice as exp (const^/X). While for a finite lattice all 
the configurations are periodic, many of them become aperiodic in the infinite volume 
limit. 

As mentioned above, except §^2 all the other sets §d contain exclusively periodic con- 
figurations. The set Spct (of plaquette-chessboard configurations) contains configurations 
built out of elementary plaquettes with filled sites, forming a square lattice according to 
the same rules as filled sites form a square lattice in the chessboard configurations from 
Scb- The remaining sets §d consist of configurations that have a quasi-one-dimensional 
structure. That is, they are built out of completely filled and completely empty lattice 
lines of given slope. Such a configuration can be specified by the slope of the filled lat- 
tice lines of the representative configuration and the succession of filled (f) and empty 
(e) consecutive lattice lines in a period. For instance, the representative configuration 
of Srfi (see Fig. E} is built out of filled lines with the slope 2 and, in the period, two 
consecutive filled lines are followed by two consecutive empty lines, which is denoted 
(2;2/, 2e). Similar description of the remaining quasi-one-dimensional configurations is 
given in Fig. El 

The coefficients {aj} for which the fourth order potentials, (^ff/'^''™°"y4) _|_ ^^^^^ 
become m-potentials are given in the tables collected in the Appendix 
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Figure 11: The {fi = 0)-phase diagram of the effective Hamiltonian ^^/^'■"^*°"^(4)_ fpj^g numbers 
in curly brackets, displayed by the symbols of open domains §d, denote the T-plaquette config- 
urations that minimize the m-potential in S_d- The numbers in curly brackets, displayed by the 
boundary-line segments or by the arrows pointing towards boundary segments (or their crossing 
points) identify the additional minimizing T-plaquette configurations. For more comments see 
the text. The boundary-line segments can be determined by means of their intersection points: 



a = (9/2,12), b = (-5/2,-2), c = (23/2,-2), d = (-5/2,-6), e = (23/2,-6), and by the 
slope 1 of the boundary of S^i, and the slope —1 of the boundary of S^i- 



Figure 12: The {fi = 0)-phase diagram of the effective Hamiltonian (ff^g'^"")*^^). The boundary- 
line segments can be determined by means of their intersection points: a = (5/2, 9), b = (1/2, 5), 
c = (9/2, 5), d = (1/2, 1), e = (9/2, 1), and by the slope 1 of the boundary of §^,1, and the slope 
— 1 of the boundary of S^i- For more explanations see the description of Fig. 111! 
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Figure 13: The (e = 0)-phase diagram of the effective Hamiltonian The 
boundary-hne segments can be determined by means of their intersection points: a = (1/2, — 1), 
b = (7/2,-3/2), c = (15/2,-3/2), and by the slope —1 of the boundaries of Si, 82, §3. For 
more comments see the description of Fig. 1111 




Figure 14: The (e = 0)-phase diagram of the effective Hamiltonian (i/gg'^°'^)(^). The boundary- 
line segments can be determined by means of their intersection points: a = (—5/2, —1/2), b = 
(11/2, -1), c = (7/2, -2), d = (3/2, -5/2), e = (15/2, -5/2), f = (19/2, -3), g = (11/2, -7/2), 
h = (7/2, —4), and by the slope —1 of the boundaries of §1 and S3. For more comments see the 
description of Fig. 
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Figure 15: Representative configurations of the sets Sd of ground-state configurations. The 
remaining configurations of can be obtained by applying the symmetries of H to the displayed 
configurations. As a representative configuration of §^2) we show a configuration with one defect 
line (the dashed line). The continuous line is a guide for the eye. For more comments see the 
text. 
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E (Tables EHIII)- 



Finally, a remark concerning ground-state configurations at the boundaries between 
the open domains is in order. Let §d and §d' be two domains of the considered 
phase diagram, sharing a boundary. At this boundary, the set of the minimizing T- 
plaquette configurations contains always the subset §td U Std', but it may contain also 
some additional T-plaquette configurations of minimal energy. Consequently, the set of 
the ground-state configurations at the boundary contains always the subset §d U S^' 
and typically a great many of other ground-state configurations, whose number grows 
indefinitely with the size of the lattice. In the considered diagrams, the only exception 
is the boundary between §seg and §>cb, where the set of the ground-state configurations 
amounts exactly to §seg U §cb- 

3.4 Discussion of phase diagrams 

We start with a few comments concerning the phase diagrams due to the fourth order 
effective interactions. In the fourth order effective interaction, the parameters uj and e 
control the strength of n.n. and n.n.n. interactions, respectively. The n.n. interaction 
is repulsive and favors the chessboard configurations if u; > 9/2, for fermions {uj > 5/2, 
for bosons). In the opposite case (attraction) it favors the S+, S- configurations. In 
turn, the n.n.n. interaction is attractive if e > 3, for fermions {e > 5, for bosons), and 
then it reinforces the tendency towards chessboard and uniform, S+, S^, configurations. 
When it becomes repulsive, it frustrates the n.n. interaction: the more negative it is the 
larger u is needed to stabilize the chessboard configurations and the smaller u is needed 
to stabilize the 5"+, configurations. Therefore, whatever the value of e is, there is a 
sufficiently large u such that {e, u) e §cb and a sufficiently small uj such that (e, u) G §seg- 
Apparently, the fourth order phase diagrams in the cases of hopping fermions and 
hopping bosons are quite similar if we compare the geometry of the phase boundaries 
and the sets of ground states. The main difference is in the domain occupying the cen- 
tral position: in the case of fermions the ground-state configurations are diagonal-stripe 
configurations, S^^, while in the case of bosons this is the set Spcb of plaquette-chessboard 
configurations. These two sets of configurations are ground-state configurations of the 
corresponding fourth-order effective interactions with the site, n.n., and n.n.n. terms 
dropped, which corresponds to the points (9/2,3) G Sdd and (5/2,5) G Spcb- This means 
that the fourth order phase diagrams can be looked upon as the result of perturbing an 
Ising-like Hamiltonian, that consists of pair interactions at distance of two lattice con- 
stants and plaquette interactions, by n.n. and n.n.n. interactions. More generally, it is 
easy to verify that the set of ground-state configurations of the Ising-like Hamiltonian, 



amount to Sdd if X > 0, and to Spcb if X < 0. These ground states are stable with 
respect to perturbations by n.n. and n.n.n. interactions, if the corresponding interaction 
constants are in a certain vicinity of zero. Otherwise, new sets of ground states, shown 
in Fig. ^2 and Fig. emerge. 

The basic question to be answered, before discussing the phase diagrams due to 
the complete interaction, refers to the relation between these phase diagrams and the 



{x,y)3 
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diagrams due to the truncated effective interactions, obtained in the previous section. 

Firstly, we note, by inspection of the phase boundaries in Fig.[TTland Fig.[T21 that the 
T-plaquette configurations {45, 58} are missing in the phase diagram corresponding to 
fermion mobile particles. In turn, the T-plaquette configurations {41,62} are missing in 
the phase diagram corresponding to boson mobile particles. Consequently, the ground- 
state configurations of Spcb cannot appear not only in the fourth order fermion phase 
diagram but also in the fermion phase diagram of the complete interaction. Similarly, 
the ground-state configurations of S^^ do not appear in boson phase diagrams. 

Secondly, by adapting the arguments presented in [2111101, we can demonstrate, see for 
instance [HI], that if the remainder R^^'^ is taken into account, then there is a (sufficiently 
small) constant to such that for t < the phase diagram looks the same as the phase 
diagram due to the effective interaction truncated at the fourth order, except some regions 
of width 0{t^), located along the boundaries between the domains, and except the domain 
§d2- For t < to and each domain §>d ^ §^2, there is a nonempty two-dimensional 
open domain Sg" that is contained in the domain and such that in S^g" the set of 
ground-state configurations coincides with 

Now, consider our systems for specified particle densities, pe = Pi = 1/2- According 
to the phase diagrams in the hole-particle symmetry plane (/i = 5 = 0), shown in Fig. ITD 
and Fig. [T21 if (£,<-i;) is in S^^, then the ground state is a phase-separated state, where 
ion configurations are mixtures of 5"+ and S- configurations, called the segregated phase. 
Another phase-separated state, where the ion configurations are mixtures of 5"+, S'_, S"^^, 
and S°j^ configurations, exists at a line that is a small (0(t^)) distortion of the boundary 
between the domains S^eg and Scb jE^I- In all the other domains §g of these diagrams 
the ground-state phase exhibits a crystalline long-range order of ions. 

By the properties of the fourth order phase diagrams, whatever the value of e is, there 
is a sufficiently large uj such that the systems are in the chessboard phase {{e,<jj) G 
and a sufficiently small oj such that the systems are in the segregated phase ((e, uS) G S^^,). 
In between, the systems undergo a series of phase transitions and visit various striped 
phases. Consider first unbiased systems = 0). Then, it follows from the phase diagram 
in Fig. ^2 that the fermion system visits necessarily the diagonal striped phase, described 
by the configurations in S^^. In turn, the phase diagram in Fig. IT^ implies that the boson 
system has to visit the dimerized-chessboard phase, described by the configurations in 
Srfi, then the vertical/horizontal striped phase, described by the configurations in §^2, and 
after that a dimerized vertical/horizontal striped phase, described by the configurations 
in §t,i. These scenarios are preserved, if e G (—2 + O(t^), 12 — 0(t^)) in the fermion case, 
and e G (— oo, 1 — 0(t^)) in the boson case. 

Note that the same sequence of transitions, as found in the unbiased boson system, 
can be realized by the fermion system if e is sufficiently small [e < —6 — 0(t^)). 

For somewhat larger values of £, £ G (—6 + 0(t^),— 2 — 0(t^)), on its way from 
the chessboard phase to the segregated phase the fermion systems visits the dimerized 
chessboard phase, then the diagonal striped phase, and after that the dimerized ver- 
tical/horizontal striped phase. The vertical/horizontal striped phase is closer to the 
segregated phase than the diagonal striped phase. In case the both kinds of stripe phases 
(vertical/horizontal and diagonal) appear in a phase diagram, such a succession of striped 
phases is perhaps generic, see also [20, ^131. 

For sufficiently large e, the tendency toward the chessboard and the uniform config- 
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urations is so strong that the both systems jump directly from the chessboard phase to 
the segregated phase. 

Above, we have described the states of the considered systems for typical values of 
the control parameter uj. Only in small intervals (whose width is of the order of 0{t^)) 
of values of cu, about the transition points of the diagrams in Fig. ^2 and Fig. ^1 the 
states of the systems remain undetermined. 

For completeness we present also ground-state phase diagrams of the both systems off 
the hole-particle symmetry plane, for e = 0, see Fig. ^[EI Away from the hole-particle 
symmetry plane new phases appear. In particular we find the phases described by the 
quasi-one-dimensional configurations §i, §2, and S3, well known from the studies of the 
phase diagram of the spinless Falicov-Kimball model [IHl- In the boson phase diagram. 
Fig. El there appears the domain §^2 that is not amenable to the kind of arguments used 
in this thesis. 
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4 Influence of nearest-neighbor anisotropy on axial 
striped phases 

In this Section^ we show that any anisotropy of nearest-neighbor hopping ehminates the 
7r/2-rotation degeneracy of the dimeric and axial-stripe phases and orients them in the 
direction of a weaker hopping. Moreover, due to the same anisotropy the obtained phase 
diagrams of fermions show a tendency to become similar to those of hardcore bosons. 



4.1 The effective interaction up to the fourth order 



We consider the model similar to (fTSj) : only a n.n. hopping is taken into account and 
we allow for its anisotropy. We set th = t, and = ^/lt, with < 7 < 1. In this 
case, up to a term independent of the ion configurations and the chemical potentials, the 
ground-state energy expansion ()43j) reads: 
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in the case of hopping hardcore bosons. As previously, the remainders, ( R^ j and 

{Rs) , are independent of the chemical potentials and the parameters W and e, but in 
this case, they depend on 7, and collect those terms of the expansions that are propor- 
tional to t^™, with m = 3, 4, . . .. The expressions (jHHI) and (jHTIjl refer to (ffH|l and (ffH|) for 
7 = 1, respectively. 

In the previous Section we have obtained the ground-state phase diagrams, according 
to the fourth-order isotropic effective Hamiltonians (7 = 1), for hole-particle symmetric 

^This section is based on [5^. 
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systems = 0) with a weak (of fourth order) n.n.n. subsidiary interaction, and for 
unsymmetrical systems without the subsidiary interaction (e = 0), see the top phase 
diagrams in Figs. ITTIIT^ Second-order phase diagrams of the isotropic case consist ex- 
clusively of phases whose configurations are invariant with respect to 7r/2-rotations, with 
macroscopic degeneracies (i.e. the number of configurations grows exponentially with the 
number of sites) at the boundaries of phase domains. Stripe phases, whose configurations 
are not invariant with respect to 7r/2- rotations, appear on perturbing the second-order 
phase diagrams by the fourth-order isotropic interactions. Here, we would like to observe 
the influence of a weak anisotropy of n.n. hopping on these stripe phases. Since we are 
working with truncated effective Hamiltonians, we have to assign an order to the devia- 
tion of the anisotropy parameter 7 from the value 1 (corresponding to the isotropic case). 
Therefore, we introduce an anisotropy order, a, and a new anisotropy parameter, jSa- 



The orders of the deviations can be neither too small, not to modify the second-order 
effective Hamiltonians, nor too high, to effect the considered effective Hamiltonian of the 
highest order. Since here, the highest order of the effective Hamiltonians is /c = 4, the 
weakest admissible deviation from the isotropic case corresponds to the highest anisotropy 
order a = 2. Then, we can consider an intermediate deviation, i.e. < a < 2. The 
strongest deviation, i.e. a = 0, is not admissible, since it modifies the second-order 
effective Hamiltonians. 

4.2 The smallest deviation from the isotropic case 

/ f\ (4) 

In the sequel, we drop the arguments of ground-state energies, that is we set I Eg ) (/i) = 



yEgj , etc. In the case of the smallest deviation from the isotropic case the fourth-order 
effective Hamiltonian for fermions reads: 



while for hardcore bosons only the first term, representing the isotropic fourth-order 
effective Hamiltonian, has to be changed properly. 

Apparently, the effective Hamiltonians of order zero and two are isotropic, and there is 
no difference between the cases of hopping fermions and hopping bosons (see Figs. I^ ITUI) . 
To observe the effect of the fourth-order anisotropy term of (j^Tj) the construction of 
the fourth-order diagrams has to be carried out again. As in the isotropic case, this is 
facilitated by introducing new variables, u, 6, and e given by (jHIl), and rewriting the 
fourth order effective Hamiltonian in the form. 
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( A 

with analogous expressions in the bosonic case, and with the isotropic potentials ( if^ j 
given by (jHH) and (jH3j) . respectively. 

7=1 

We would like to get an idea of the phase diagram in the space of the four energy 
parameters (tu, e, ^2-, ^\ that appear in the Hamiltonian. Due to the fact that the domains 
occupied by the phases are polyhedral sets, this goal can be achieved by studying phase 
diagrams in two-dimensional hyperplanes. Of particular interest are those hyperplanes 
that result from intersecting the four-dimensional space by hyperplanes 5 = and (i^ = 
const, and by hyperplanes e = and /?2 = const. A collection of such sections for suitable 
values of /32, forming a finite increasing from zero sequence, enable us to observe how the 
anisotropy effects our system if it is hole-particle symmetric and if it is not, respectively. 

Specifically, in Fig. we show phase diagrams in the hole-particle-symmetric case, 
while in the absence of the hole-particle symmetry, the phase diagrams are shown in 
Fig- El ^ for fermions, and in Fig. for bosons. 

The phases that appear in these phase diagrams can conveniently be described in 
terms of the phases found in ground-state phase diagrams of the isotropic systems, with 
Hamiltonian H given by (fTHj) (see Fig. [T3j). The phases of anisotropic systems, considered 
here, either coincide with or are simple modifications of the isotropic phases. In fact, in 
Fig.Hnionly representative configurations of the phases of isotropic systems are displayed. 
The remaining configurations can be obtained by means of the spatial symmetry oper- 
ations of isotropic Hamiltonian, like translations and rotations by n/2. The numbers in 
curly brackets, placed by the symbols of phases, stand for the numbers of the T-plaquette 
configurations (according to Fig. IHUj) that are obtained by restricting the configurations 
of a phase to a T-plaquette. By the same symbols as the phases we denote also their 
domains. 

Among the phases of isotropic systems, we can distinguish a class of dimeric phases, 
Srfi, . . . , Srf4, and a class of axial- stripe phases, S„i, §„2, and 8^,3. The configurations of 
dimeric phases consist of isolated pairs of n.n. occupied sites (in the sequel called the 
(timers). In the configurations of axial-stripe phases, the ions fill completely some, parallel 
to one of the axes, lattice lines, so that a periodic pattern of stripes is formed. Out of 
the dimeric or axial-stripe configurations of Fig. only those with dimers or stripes, 
respectively, oriented vertically appear in the phase diagrams of anisotropic systems. 
Such a restricted phases are marked in anisotropic phase diagrams by the additional 
superscript, v. 

In comparison with the isotropic phase, a new domain §^4 appears in the phase 
diagrams shown in Fig. El Similarly to §^2, its degeneracy grows indefinitely with the 
size of the lattice. This phase consists of three classes of periodic configurations of dimers. 
In the class Sd4a, the elementary cell can be chosen as a parallelogram whose two sides 
of length 3 are parallel to dimers (which are vertical or horizontal). If the dimers are 
oriented vertically, then the other two sides have the slope 1/2 and the length ^/E. By 
reflecting an elementary cell of §^40 in a lattice line passing through its side that is 
parallel to dimers, we obtain an elementary cell of the class §^4;,. In the third class, SdAc, 
an elementary cell can be chosen as a rhomb formed by the centers of dimers, with the 
sides of length vTO. Two configurations, one from Sdia and one from 8^4;,, having the 
same kind of dimers (vertical or horizontal), can be merged together along a "defect line" 
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Figure 16: The case of the weakest anisotropy and of hole- particle symmetric systems (fj, = 0). Phase diagrams of 
(^^t) (given by 1931 ^ — left column, and — right column, for an increasing sequence of values of (^2- The 

representative ion configurations of the displayed phases are shown in Fig. 1151 f for more comments see text). For fermions, 
the critical value is /32 = 7, while for bosons it is /32 = 2. The equations defining the boundary lines of the phase domains 
are given in Tab. [Tl and Tab. 1^ of Appendix D, while the corresponding zero-potential coefHcients {cti} — in Tab. I12l — 
Tab. ll9l of Appendix E. In the bottom diagrams, the regions surrounded by dotted circles, are reconsidered in the case of 
an intermediate anisotropy. 
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Figure 17: The case of the weakest anisotropy, off the hole-particle symmetry, with e = 0. The phase diagram of 
(^^t) (given by for an increasing sequence of values of 132 ■ The representative ion configurations of the displayed 

phases are shown in Fig. 1151 ffor more comments see text). The critical values are: (32 = 1, /32 = 2, and (32 = 3. In 
the blank region of the phase diagram for the critical (32 = 3, the following T-plaquette configurations have the minimal 
energy: 26,52,53,68,81,88,101,116,117,143 (see Fig. 1301 . By means of these T-plaquette configurations one can construct 
1 Sddi and many other configurations. The equations defining the boundary lines of the phase domains are given in 
Tab. Islof Appendix D, while the corresponding zero-potential coefficients {oi} — in Tab. 1201 - Tab. l26l of Appendix E. In 
the bottom diagram, the regions surrounded by dotted ellipses are reconsidered in the case of an intermediate anisotropy. 
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Figure 18: The case of the weakest anisotropy, off the hole-particle symmetry, with e = 0. The phase diagram of 

(-ffy)'*' (given by I93i 'l: the isotropic diagram (/32 = 0) and the anisotropic diagram (/32 = 3). The representative ion 
configurations of the displayed phases are shown in Fig. 1151 (for more comments see text). Here, no critical values have 
been detected. The equations defining the boundary lines of the phase domains are given in Tab. HI of Appendix D, while 
the corresponding zero-potential coefficients — in Tab. I27l and Tab. I28l of Appendix E. In the right diagram, the 

regions surrounded by dotted ellipses are reconsidered in the case of an intermediate anisotropy. 

parallel to dimers (dashed line in Fig. ITTH) without increasing the energy. In this way, a 
numerous family of configurations can be constructed, with the number of configurations 
growing like exp {const\^). 




Fi gUre 19: The representative configurations of the set SJJ^ whose degeneracy grows like exp (const v'A). The left 
configuration is an example of configurations with defect lines. Here, there is one defect line (the dashed line, the continuous 
line is a guide for the eye): the vertical lattice line separating a periodic configuration of vertical dimers from its vertical 
translate by one lattice constant. The right configuration is a periodic configuration of dimers. For more comments see 
text. 

It follows from the polyhedral shape of the phase domains, that the set of values of /?2 is 
partitioned into open intervals, where the boundary lines of phase domains do not change 
their direction, only their distance to the origin varies in an affine way. The boundary 
points of these open intervals are special values of the anisotropy for the phase diagrams, 
and in the sequel we call them the critical values. As a critical value of anisotropy is 
approached, some boundary lines merge into a line or a point, which results in disappear- 
ance of some phase domains. And vice versa, some points and lines break off, creating 
new phase domains. 

In particular, in the hole-particle symmetric case it can be inferred from (Fig. ITT)|l that 
at least up to /?2 = 10 - for hopping fermions, and at least up to = 5 - for hardcore 



48 



bosons, there is only one critical value of /?2- For fermions, it amounts to (32 = 7, where 
the phase S^^ disappears, while for hardcore bosons it is /?2 = 2, where the phase Spcb 
disappears. 

Off the hole-particle symmetry, for bosons, there are no critical values of (32-, at least 
up to (32 = 3. On the other hand, for fermions and for (32 < Q there are as many as three 
critical values of ^2- In increasing order, the first is P2 = 1, where the phases S^i, SJJ^, 
and §^2 appear. The second is P2 = 2, where the phases SJ^g, §^^4, and S^g appear. And 
the last one is /32 = 3, where S^^ is replaced by §^3- 

If there are no more critical values of the anisotropy parameter (32, the shape of 
the fourth-order phase diagrams for any value (32 larger than the greatest considered in 
Figs. El ^1 remains the same. On increasing the anisotropy parameter, the diagrams 
undergo only some translations. To verify whether for the values of (^2 larger than those 
considered in this section new critical values do appear, we proceed to investigating 
stronger, i.e. intermediate, deviations from the isotropic case. 

4.3 The intermediate deviation from the isotropic case 

Now, the fourth-order effective Hamiltonian for fermions assumes the form: 

/ \ (4) / \ (4) t^^'^ 

and for hardcore bosons an analogous formula holds true. Obviously, the phase diagrams 
in the zeroth and second orders remain the same as in the isotropic case, described 
above. Therefore, we proceed to constructing the phase diagram in next order, which is 
(2 + a)-order with < a < 2, and the corresponding effective Hamiltonian reads: 

^(2+a)^^(2)_^^t^ ^ S^Sy, (95) 

(2) 

where Eg stands for the, common for fermions and bosons, second-order effective Hamil- 
tonian. For the reasons given in the previous subsection, we consider a neighborhood of 
the point W = — 2t^, /i = 0, e = 0, where the energies of all the configurations are equal. 
In this neighborhood it is convenient to introduce new variables, S', e', and u', 

^ = e+''5\ e = t^+^e', W = -2f + t^+'^u', (96) 

and rewrite the expansion up to the order 2 + a: 

J ^1 




{x,y)i {x,y)2 {^,y)i,v 



(97) 



where 



^ {x,y)i {x,y)2 {x,y)i,v 
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The summations in the primed sums run over a plaquette P2. The plaquette potentials 
j^g^yg minimized over all the plaquette configurations. As in the previous 

case, we are interested in phase diagrams of hole-particle symmetric systems {6' = 0) 
or unsymmetrical systems with e' = 0. It turns out that for such energy parameters 
and plaquette configurations the potentials Hp^"'^ are m-potentials. The resulting phase 
diagrams are shown in Fig. EOland Fig. ED We note that each of the points a: to' = Pa, 




Figure 20: The case of an intermediate anisotropy 
(0 < a < 2) and of hole-particle symmetric systems {fj, = 
0). The phase diagram (common for hopping fermion and 
hardcore boson systems) of Hp^'^K The coordinates of 
point a are uj' = Pa, e' = Pa- The boundary lines of S^ji 
from left to right, are: e' = a;' and e' = —w' + 2/3a- The 
boundary between Sseg and S^t is = /3a- 



Figure 21: The case of an intermediate anisotropy 
(0 < a < 2), off the hole-particle symmetry, with e' = 0. 
The phase diagram (common for hopping fermion and 
hardcore boson systems) of Hp^^"'. In the (cj', 5')-plane, 
b= (0,0), c+ = (2/3a,/3a), c- = (2/3a,-/3a). The bound- 
ary lines of S_, from left to right, are: <5' = 0, 5' = — ^ 
and <5' = —uj'+l3a. The boundary lines of §+ are obtained 
by changing S' — <5'. 



e' = (3a, b: u' = 0, 6' = 0, c": u' = 2(3a, 6' = ~f3a, and c+: u' = 2(3a, 6' = (3a, is 
the coexistence point of three periodic phases. That is, the only plaquette configurations 
minimizing the potential Hp^""^ at such a point are those obtained by restricting the 
configurations of coexisting phases to a plaquette (these plaquette configurations are 
shown in Fig. U^i 

Now, following our recursive procedure of constructing phase diagrams to some order, 
we are ready to investigate the effect of fourth-order interactions. As in earlier steps, it is 
enough to consider neighborhoods of the coexistence points a, b, and c~ of Figs. 1201 1^ 
(the diagram in a neighborhood of c+ can be obtained from that around c~ by a symmetry 
operation). 

A. A neighborhood of a 

Here 6' = 0, hence the systems are hole-particle invariant. A convenient change of 
variables is: 

uj' = Pa + e' = (3a + t'-'^e. (99) 



50 




Figure 22: The case of an intermediate anisotropy 

(0 < a < 2) and of hole-particle symmetric systems 

/ f \ (4) 

{fi = 0). The phase diagram of [E-'g] in a neighbor- 
hood of point a (w' = f3a, e' = /3a) (see formula llOUi '). 
The equations defining the boundary lines of the phase 
domains are given in Tab. ISlof Appendix D. 



Figure 23: The case of an intermediate anisotropy 
(0 < a < 2) and of hole-particle symmetric systems 
{/i = 0). The phase diagram of (Bg)'*' in a neighbor- 
hood of point a {uj' = f3a, e' = /3a) (see formula llOUi '). 
The equations defining the boundary lines of the phase 
domains are given in Tab. 151 of Appendix D. 



In these variables, the fourth-order effective Hamiltonian reads: 



El) 



(4) 



^2+a 



P2 



(2+a) 
P2 



e'=uj'=l3a T 



(4) 



(100) 



7=1 



with a similar formula for hard-core bosons, where 

(3a 



H 



{2+a) 
P2 



5'=0 




(101) 



In the previous order, it has been established that the only plaquette configurations 

are the ones obtained by restricting to a plaquette P2 the pe- 



mmimizmg Hp^ 



riodic configurations §+, §c6, and §^2- Let us denote this set of plaquette configura- 

/ r\ (4) 

tions by Sp^ • Consequently, the minimization of the fourth-order potentials ( Hj^ j 



6\(4) 



7=1 



7=1 



should be carried out only over the set of those T-plaquette configurations 



whose restriction to a plaquette P2 belongs to Sp^ . The configurations of the set Sf. are dis- 



played in Fig. EH It appears that on the set Sf., the potentials 



(4) 



M(4) 



7=1 



7=1 



are m-potentials. The obtained phase diagrams, independent of the anisotropy parameter 
/3a, are shown in Fig. |221aiid Fig. 1221 



B. A neighborhood of b 

Here the system is not hole-particle symmetric and e' 
variables is: 



6' = r-^d. 



0. A convenient change of 



(102) 
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Then, the fourth-order effective Hamiltonian takes the form, 



E 



(4) t 



2+a 



P2 



4) 



(4) 



7=1 



where 



H 



(2+a) 



P2 



5'=0 
e'=uj'=0 



4e' 



SxSy 



{x,y)i 



with a similar formula for hard-core bosons. The minimum of Hp^""^ 



S'=0 



(103) 



(104) 



is attained 



at the configurations belonging to Sp^, i.e. the plaquette configurations obtained by 
restricting the periodic configurations §+, S_, and a plaquette P2. Let be the 

corresponding set of T-plaquette configurations (there are no vertical pairs of n.n. sites 

A (4) 



occupied by one ion). This set is shown in Fig. |nH Here the potentials (^-^t) 
not the m-potentials. The obtained phase diagrams are shown in Figs. EH ESI 



are 



7=1 



C. A neighborhood of c~ 

Here the system is not hole-particle symmetric and e' = 0. A convenient change of 
variables is: 



The fourth-order effective Hamiltonian reads: 



(4) t2+a 



P2 



e'=0 

5' = -l3a,Ul'=2l3a 



(4) 



7=1 



with a similar formula for hard-core bosons, where 



rr(2+a) 
^P2 



e'=0 

5' = -l3a,Lu' = 2l3a 



Pa 

4 



^ ^ i,Sx ~^ ~^ SxSy -\- 1^ 

{^,y)i,h 



The potentials Hp_^^^ 



(105) 



(106) 



(107) 



e'=0 

5' = -Pa,i0' = 2f3a 



are minimized by restrictions to a plaquette P2 of 



periodic configurations §_, §cb, and §^35 ^^at constitute the set Sp^. The corresponding 
set §^ of T-plaquette configurations consists of configurations where no horizontal pair 



of n.n. sites is occupied by two ions (see Fig. 1221) • Here the potentials [Hr 



(4) 



are not 



7=1 



the m-potentials. The corresponding phase diagrams are shown in Fig. ESI and Fig. EZI 



4.4 Discussion of the phase diagrams and conclusions 

In previous subsection, we have obtained phase diagrams according to the truncated 
effective Hamiltonians of fourth order, with an intermediate anisotropy of hopping. These 
diagrams have been constructed to see what happens to the fourth-order phase diagrams 
with the weakest anisotropy if the anisotropy parameter P2 grows beyond the values taken 
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5 



5 



Figure 24: The case of an intermediate anisotropy 
(0 < a < 2), off the hole-particle symmetry, with e' = 0. 
The phase diagram of (^5) ' ' in a neighborhood of 
point b (iv' = 0, <5' = 0), see formula OSSJ- The equa- 
tions defining the boundary lines of the phase domains 
are given in Tab. |H] of Appendix D, while the corre- 
sponding zero-potential coefficients {a^} in Tab. of 
Appendix E. 



Figure 25: The case of an intermediate anisotropy 
(0 < a < 2), off the hole-particle symmetry, with e' = 0. 
The phase diagram of [Eg) in a neighborhood of point 
b (u)' = 0, 5' = 0)), see formula llU3i . The equations 
defining the boundary lines of the phase domains are 
given in Tab. ISlof Appendix D, while the corresponding 
zero-potential coefficients {oi} in Tab. lijOl of Appendix 
E. 




Figure 26: The case of an intermediate anisotropy 

(0 < a < 2), off the hole-particle symmetry, with e' = 
( f\ 

0. The phase diagram of ( Ei, j in a neighborhood 
of point {w' = 2/3a, S' = —/3a), see formula llU6i . 
The equations defining the boundary lines of the phase 
domains arc given in Tab. |7|of Appendix D, while the 
corresponding zero-potential coefficients {a;} in Tab. l3li 
of Appendix E. 



Figure 27: The case of an intermediate anisotropy 
(0 < a < 2), off the hole-particle symmetry, with e' = 
0. The phase diagram of [Eg) in a neighborhood 
of point {uj' = 2/3a, (5' = —Pa), see formula 11061 . 
The equations defining the boundary lines of the phase 
domains are given in Tab. |7|of Appendix D, while the 
corresponding zero-potential coefficients {oi} in Tab. l32i 
of Appendix E. 
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into account in subsection 4.2. Let us consider firstly the diagrams in a neighborhood of 
point a. Apparently, up to a rescaling the phase diagram in Fig. |221 looks the same as 
the part in the dotted circle of the fermionic phase diagram in Fig. ^1 Similarly, up to a 
rescaling the phase diagram in Fig. 1221 looks the same as the the part in the dotted circle 
of the bosonic phase diagram in Fig. ^1 The same remarks apply to the phase diagrams 
in neighborhoods of points b and c~. The phase diagram in Fig. |21] reproduces the part 
in the upper dotted ellipse in Fig. El while the diagram in Fig.|2Sl^ the part in the upper 
dotted ellipse in Fig. Then, the phase diagram in Fig. 123 reproduces the part in the 
lower dotted ellipse in Fig. El while the diagram in Fig. |2Z1- the part in the lower dotted 
ellipse in Fig. El Therefore, we conclude that in the case of the weakest anisotropy 
in the fourth-order effective Hamiltonians, we have determined all the critical values of 
the anisotropy parameter ^2- For the both kinds of hopping particles, if P2 exceeds the 
greatest critical value, the obtained phase diagram undergoes only a translation (varying 
with ^2)- 

Now, the basic question to be answered is concerned with the relation between these 
fourth-order phase diagrams and the phase diagrams of quantum systems described by 
Hamiltonians H, given by (f7^ . By adapting the arguments presented in Refs. jSHllinj, we 

can demonstrate, see for instance Ref . [M| , that if the remainders, (^Rs^ and , 
are taken into account, then there is a sufficiently small to such that for t < to the phase 
diagrams of quantum systems look the same as the phase diagrams according to the 
effective Hamiltonians truncated at the fourth order, except some narrow regions, of width 
0{t^) (at the dia grams displayed above), located along the phase-domains boundaries, 
and except the domains §^2 §^4. For t < to and for each domain $d, which is different 
from a'^d §^4, there is a nonempty two-dimensional open domain S'g' that is contained 
in the domain and such that in the set of ground-state configurations coincides 
with Moreover, in comparison with the critical values of the anisotropy parameter 
f32, determined according to the fourth-order effective Hamiltonians, the corresponding 
critical values of the quantum systems described by Hamiltonians H differ by 0{t^), i.e. 
for a quantum system 7 = 1 — /?2i^^ + 0(t^). 

Remarkably, in the fourth-order the hole-particle symmetric phase diagrams of fermions 
and of hardcore bosons are geometrically similar. That is, a phase diagram of hard core 
bosons, with any > 0, can be obtained from a phase diagram of fermions, with ^ 5, 
by the translation whose vector reads: u = —7, e = —3, and (32 = —5. The existence of 
this translation vector is related to the fact that for both kinds of systems there is one 
critical value of the anisotropy parameter. Additionally, for fermions with /?2 < 7, it is 
necessary to replace the phase §dd in the central domain by Spcfe- Off the hole-particle 
symmetry, the relation between the bosonic and fermionic phase diagrams is not that 
close. For fermions, there are three critical values of f32, while for hardcore bosons there 
is no critical values. Thus, the system with hopping hardcore bosons is less sensitive 
to the anisotropy of hopping, than the system with hopping fermions. Nevertheless, if 
e = 0, then the phase diagrams of both kinds of systems are topologically similar, except 
that in the bosonic phase diagrams the phases S^g and §2 are missing. However, we know 
from Ref. [HI] that this deficiency can be removed by switching on the n.n.n. interactions 
with negative e. 

In the fourth-order effective Hamiltonians (jH^, the weakest anisotropy of n.n. hopping 
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assumes the form of a fourth-order attractive n.n. interaction in vertical direction (i.e. 

the direction of a weaker hopping). This interaction favors n.n pairs of occupied or empty 
sites that are oriented vertically. As a result, the dimeric and axial-stripe phases oriented 
vertically are stabilized for any value of the anisotropy parameter P2, while Spcb and §m 
are replaced by §^2 above a critical value of /32- Note however, that at any higher order 
2k, k — 3,4... the weakest anisotropy of n.n. hopping will cause the same effects, in 
the effective Hamiltonians as well as in the corresponding phase diagrams. This implies 
that in the quantum systems described by H, arbitrary small anisotropy of n.n. hopping 
orients the dimeric and axial-stripe phases in the direction of a weaker hopping. 
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5 Influence of next-nearest-neighbor anisotropy on 
diagonally-striped phases 



In the previous Section the influence of n.n. anisotropy of hopping intensity on axial 
striped phases, was investigated. At the regime where stripes are stable, we have proved 
rigorously that for both systems, of hopping fermions and hard-core bosons, an arbitrarily 
small anisotropy of n.n. hopping orients the axial striped phases in the direction of a 
weaker hopping. The analogous, arising naturally, question is how the anisotropy of 
the next-nearest-neighbor (n.n.n.) hopping influences the degeneracy of diagonal-striped 
phases. To answer this question, we investigate the influence of n.n.n. hopping on the 
phase §dd (Fig. HH), whose stability was proved for fermions in Section 3. 

In this Section^ we introduce a next-nearest-neighbor hopping, small enough not to 
destroy the striped structure and examine rigorously how the presence of the next-nearest- 
neighbor hopping anisotropy reduces the 7r/2-rotation degeneracy of the diagonal-striped 
phase. The effect appears to be similar to that in the case of anisotropy of the nearest- 
neighbor hopping: the stripes are oriented in the direction of the weaker next-nearest- 
neighbor hopping. 



5.1 The effective interaction up to the fourth order 



We consider our model in the case when only next-nearest-neighbor hopping is anisotropic, 
and nearest-neighbor hopping amplitude remains the same in different directions, i.e. 
th = tv = t. For simplicity 6 = 0. In this case, the ground-state energy expansion 
with the expansion terms up to the fourth order, i.e. a + 6 + c ^ 4, reads: 



Esifi) =Ef + 



+ 



+ 



+ 



+ 



4 16 16 ^ + + 8 
4+16 8 ^ + + ^ 16 + 8 + - 

1^2 + _ 3^2 /2t2 + t t_) _ _ 3^2 ,2 

4" 16 8^" ^ ' 16" 8+" 



{x,y)i 
{^,y)2- 



1 4 1 2 
-t^ - -t^t+t. 

8 8 + 



16 + " 



E 3 2 2 X ^ 3 2 2 X ^ 



{^,y)3 



{x,y)A,+ 



{x,y)4 



1 A ^ ^ 1 A ^ ^ 3 2 \ 3 2 \ ^ 

+ fi^+ 2^ ^^^y + s^- ^^^y + fi^ V + ^- 2^ 



{^,y)5,+ 



{^,y)5.- 



Pi 



16 



P2 



Pi 



Pi 



^This section is based on [HU], 
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(108) 



The remainder i?*^^) is independent of the chemical potentials and W ^ and collects all 
the terms proportional to with a + 6 + c = 5,6, . . .. The above expansion is 

absolutely convergent for sufficiently small t, t+ and t_, uniformly in A. In the special 
case of = t_ = t' and W = Q ii was obtained in |64j . 

In Section 3, we have obtained the phase diagram of the isotropic model without 
n.n.n. hopping, i.e. for t+ = t_ = 0. Here, our aim is to determine the influence of the 
n.n.n.-hopping anisotropy on the diagonal-striped phase $dd (see Fig.lTTj). For this job, the 
value of the n.n.n.-hopping intensities, t+, t_, cannot be too large, in order to preserve the 
phase diagram up to 4th order. On the other hand, the n.n.n.-hopping intensities cannot 
be too small, in order to appear in the fourth-order effective Hamiltonian. In an attempt 
to satisfy the both requirements, we choose the smallest n.n.n.-hopping intensities t+, 
t_, i.e. such that they do not appear in the expansion terms of order smaller than four: 



= a^t and t_ 



a_t . In this case, the effective Hamiltonian assumes the form: 



Ef (^) 



1 3 4 , 



1 2 9 4 

-t^ + — 

4 16 



+ 



4 + 16 



SxSy 



16 



E 

+ ^-^y + ^^'^+ + ^p, + Y 

p+ p- 



{x,y)3 



(109) 



P2 



i.e. the second requirement is satisfied. To answer the question concerning the infiuence 
of anisotropy of n.n.n. hopping on the degeneracy of the phase S^^, there is no need 
to consider the whole phase diagram. For t+ = t_ = 0, we fix a point, well inside the 



domain of the diagonal-striped phase Sdd, say yU = and W = —2t + 9/2t , i.e. 



UJ 



9/2 



(see Fig. Ill|). Then, with the fixed point in (/x, iy)-plain, we introduce a n.n.n. hopping 
which does not change the ground-state configurations. Calculations show that S^^ has 
the minimal energy for a = |a+| = |a_|, where —1/4 ^ a ^ 1/4 (we suppose that the 
difference between a+ and a_ is not large, so they are of the same sign). Therefore, 
with our choice of n.n.n.-hopping intensities, the first of the above two requirements can 
also be satisfied. Eventually, we fix the values of n.n.n.-hopping intensities: a+ = 1/8, 
a_ = 7a+, with 7 varying about 1 (say, ^ 7 ^ 2). 

Since all the energy parameters, except the parameter 7 of n.n.n.-hopping anisotropy, 
have been fixed, the effective Hamiltonian (jl09p depends only on 7. In the following 
subsection, we examine how n.n.n.-hopping anisotropy, 7 7^ 1, influences the degeneracy 
of the diagonal-striped phase S^^- 



5.2 Diagonally-striped phase versus n.n.n.-hopping anisotropy 

For technical reasons (see paragraph 2.4.5), it is convenient to deal with such energies of 
conflgurations that are affine functions of the parameters of the effective Hamiltonian. 
However, the effective Hamiltonian ()109|) contains the terms proportional to 7 and 7^. To 
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Figure 28: The phase diagram of £''^^^((^1,^2)5 with ^ ^ 2 and ^ ^2 ^ 4. The 
rectangle breaks down into two regions: above ^2 = 1 the configurations S^^ are the 
ground-state configurations, while are replaced by below d2 = I. The dashed line 
represents the condition d2 = d\. The representative configurations (up to translations) 
of §^^, S^^ are shown on the right. 



get rid of nonlinearities, we replace 7 and 7^ by two independent parameters di and ^2, 
respectively, with di, c?2 varying in the rectangle < < 2 and < (i2 < 4, in which the 
Hamiltonian is affine. After constructing the phase diagram in (rfi, (i2)-plane, we restrict 
it to the d2 = dl curve. 

To compare the energies of configurations, we rewrite E'g\di,d2) as the sum. 



Ef\d„d2) = '-J2H?, (110) 

T 

over T-plaquettes. The potential h!^^ is of the form: 

{^,y)2,+ {x,y)2- {x,y)3 

3 v^" 3 , 5 

+ 1^ L + I^^i L 'P.- + ^ L (111) 
p+ p- P2 

where S5 is the central site of a T-plaquette (Fig. [7j). 

Unfortunately, the potential (jlllj) is not an m-potential in the rectangle of consid- 
ered values of di and d2. Therefore, we introduce appropriate zero-potentials given in 
Appendix A. In order to obtain the phase diagram, we have to compare the energies of all 
the possible T-plaquette configurations. The zero-potential coefficients a needed for this 
are given in Tab. 1331 in the Appendix. We provide their values only at certain generating 
points, since we can assume that the coefficients a are affine functions of parameters 
(^1,^2). 

The phase diagram of E^^^ is shown in Fig. |2H1 The rectangle of considered points 
{dl, (^2) breaks down into two domains. In the lower one, where 7 < 1 and t_ < t+, it is 
the phase S^^, with stripes being parallel to the direction of t_-hopping, that is stable. 
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The analogous situation is in the upper domain, where 7 > 1 and t_ > t+: the stable 
phase, S^^, consists of stripes oriented along t+-hopping. At 7 = 1, we have the isotropic 
phase §dd whose configurations consist of S^^ and S^^. 

So we see that, at least for the truncated effective Hamiltonian, switching on of a 
n.n.n.-hopping anisotropy reduces the rotational degeneracy of diagonal-striped phases: 
they become oriented in the direction of the weaker hopping. 

This result is similar to that described in Section 4, where the influence of n.n.- 
hopping anisotropy on axial-striped phases was investigated. In that case, not only for 
a truncated effective Hamiltonian but also for the corresponding quantum one, it was 
proved that for any nonzero value of n.n. -hopping anisotropy the rotational degeneracy 
of axial-striped phases is reduced by making them oriented along the direction of the 
weaker n.n. hopping. Now in turn, the natural question is whether the conclusions 
we arrived at, concerning n.n.n.-hopping anisotropy, hold true for the quantum model, 
described by Hamiltonian (fTHj) . Applying the arguments presented in jlOlinn], it could be 
demonstrated, that the stable phases of the obtained above phase diagram remain stable 
for the model (fTSj) . but in some smaller domains. That is, if the remainder i?^^) is taken 
into account, then there exist such a small to, that for t < the phase diagram looks the 
same for the quantum model, excepting of some narrow regions (of width 0{t) in the scale 
of the fourth-order phase diagram shown in Fig. OHj) . located along the phase-boundary 
lines. In our case that means that the breaking of the rotational symmetry occurs for 
7 = 1 + 0{t), when the n.n.n.-hopping intensities are O(t^). Unfortunately, we cannot 
claim that any non-zero n.n.n.-hopping anisotropy reduces the rotational degeneracy of 
the quantum model, as it was the case for n.n. hopping (see Section 4). Here it seems, 
at least for small n.n.n.-hopping intensities, that there is certain critical value of I7 — 1|, 
above which the degeneracy of phase §dd is reduced. 
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6 Conclusions 



In Sections 3 and 4, we have succeeded in constructing ground-state phase diagrams of 
some classical (i.e. diagonal in some basis) Hamiltonians, in the whole range of parame- 
ters (coupling constants, chemical potentials). These Hamiltonians have been derived as 
effective Hamiltonians in certain orders of quantum Hamiltonians (termed the extended 
Falicov-Kimball Hamiltonians), in the strong-coupling regime. In other words, we have 
proved that the considered classical Hamiltonians can be expressed by m-potentials, for 
arbitrary value of parameters. Such Hamiltonians are sometimes referred to as frus- 
trationless. It should be emphasized however, that, a priori, nothing guarantees that 
the method we used here will "produce" any phase diagram: one can fail to construct 
m-potentials or they might not exist. 

Those ground-state phase diagrams of effective Hamiltonians constitute the basic 
input of the quantum extension of the Pirogov-Sinai theory of phase diagrams, developed 
in j2ZI. This theory enables us to claim that not only the periodic ground states of effective 
classical Hamiltonians are also ground states of the corresponding quantum Hamiltonians, 
but that, at sufficiently low temperatures, to each periodic ground-state phase there 
corresponds a thermodynamic phase characterized by the same kind of long-range order. 
The above mentioned properties of periodic ground-states of effective Hamiltonians are 
referred to as stability with respect to higher-order effective interactions and to small 
thermal fluctuations. 

Several basic problems concerning striped phases in a system of strongly correlated 
electrons were considered. First of all, we addressed the problem of the stability of 
charge-stripe phases versus the phase separated state — the segregated phase, in strongly- 
interacting systems of fermions or hardcore bosons. There are numerous works devoted 
to this problem, where it is studied, by approximate methods, in the framework of mod- 
els relevant for experiments. Unfortunately, due to tiny energy differences involved, it is 
difficult to settle this problem by means of approximate methods, which bias the calcu- 
lated energies with hardly controllable errors of various nature. We have studied simple 
models, that by many physicists can be considered less realistic, but which, in return, are 
amenable to a rigorous analysis. Our models are described by extended Falicov-Kimball 
Hamiltonians, with hopping particles being spinless fermions or hardcore bosons. The 
ground-state phase diagrams of these models have been constructed rigorously in the 
regime of strong coupling and half-filling. In the both cases, of fermions and hardcore 
bosons, we have found transitions from a crystalline chessboard phase to the segregated 
phase via striped phases. 

After the stability of striped phases had been proven, both system (with the hopping 
particles being either spinless fermions or hardcore bosons) have been studied in order 
to show rigorously the influence of nearest-neighbor and next-nearest-neighbor hopping 
anisotropy on phase diagrams. In the case of nearest-neighbor hopping only, two main 
conclusions have been drawn. Firstly, arbitrary small anisotropy of nearest-neighbor 
hopping orients the dimeric and axial-stripe phases in the direction of a weaker hop- 
ping. Secondly, even a weak anisotropy of hopping reveals a tendency of fermionic phase 
diagrams to become similar to the bosonic ones. 

In the case of next-nearest-neighbor hopping we have shown that a weak anisotropy 
of the next-nearest-neighbor hopping reduces the degeneracy of a diagonal-striped phase. 
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it orients the stripes in the direction of the weaker next-nearest-neighbor hopping. 
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Appendix A 

Here we present T- block zero-potentials in the isotropic case (considered in Section 3), 
which are invariant with respect to 7r/2- rotations and reflections. 

k!:p (s) =Si + S3 + Sy + Sg - As^, 

k?\s) ^s, + s, + s, + Ss-As,; 

- 2S2S5 - 2S5S6 - 2S5S8 - 2S4S5, 

k!^\s) =siS5 + S3S5 + S5S9 + S5S7 - S2S4 - S4S8 - sgse - S2S6, 

kr [s) =siS3 + S3S9 + .s"7,s"9 + S1S7 - 2S4SQ - 2s2S%; 

kP{s) =SiS2Ss + S3S6S9 + SjSsSg + S1S4S7 - 2S4S5S6 " 25255^8, 

k)^'[s) =818285 + 848587 + S5S8S9 + S3S5S6 + 828385 + S1S4S5 + 858788 + S5S6S9 

- 2818284 — 284878s — 2868s8g — 2828386, 

kr {s) =818285 + 848587 + S5S8S9 + S3S5S6 + 828385 + 818485 + S58788 + 858689 

- 2828485 — 284S58S — 285S6SS — 2828586, 

k!^\8) =818286 + 828487 + S4S8S9 + 83868s + 86878s + S2S6S9 + S28384 + 81848s 

- 83S485 - 81858s - 858687 - S28589 - 848589 - 83S58S - 818586 - 828587, 

k^°\s) ^818287 + 848789 + 838S89 + 818386 + 82S389 + 818384 + 81878s + 868789 

- 81828s - 848687 - 828S89 - 838486 - 82838s - 818486 - 82878s - 848689, 
k^^\s) ^818385 + 818587 + 858789 + S38589 - 82848s - 84868s - 82868s - 828486] 

(s) =81828487 + S4878S89 + 83868S89 + 81828386 + S2838689 + 8182S384 

+ 8184878s + SgSySgSg — 8182858s — 84S58687 — 8285SS89 — 83848586 

- S2S385SS - 81848586 - 828587SS - 84858689, 

kT^\8) =81828385 + 81848587 + 85878SS9 + 838586S9 - 82S4858S - S48586SS 

- 828586SS - S2S4S5S6, 

k^'^\8) =81828488 + 5456-57-58 + -S2-S6-S8-S9 + -S2-S3-S4-S6 + 8283868s + -S1S2-S4-S6 

+ 8284878s + 84868S89 — 81828587 — 84858789 — 535558-89 ~ 51535556 

- 52535559 - 51535455 - 51555758 - 55555759; 

kT^\8) =8486878889 + 5253565859 + 5i52535456 + 5i52545758 - 5455555759 
~ -52-^3-^5-58-^9 ^ -^153545556 — 5i52555758, 

A;^^'*(s) =5355555859 + 5154555753 + 5152535555 + 5152545557 + 5455575859 

+ 5555575359 + 5253555559 + 5152535455 - 5455555359 - 5253555553 

- 5152545555 — 5254555753 — 5253545555 — 5i52545558 — 5455555753 

- 5255555359. 

In the case of n.n. hopping anisotropy (considered in Section 4), the zero-potentials 



68 



invariant with respect to 7r-rotations and reflections are of the form: 

4S5, 
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,(3: 
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,(5: 
,(6: 

Vr-j-1 
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brj-i 
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l(16 

ri/rj-i 
krp 
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A/rj-i 
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fbrj-i 

(20 



k. 
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(21 
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(22 
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(23 
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1,(24: 

l(25 

(26 
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(27' 
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(28 



k. 



(29) 



2S4S5 - 2S5S6, 
2S2S5 - 2S5S8, 

S2S4 - S4SS - ssse - S2S6, 



S2S4S5 
S2S4S5 



55=^6=^8 ~ S2S5S6 



S4S5S8) 
545558) 
545558, 



51 + S3 + S7 + Sg 

52 + 58 - 2S5, 
54 + 56 - 2s5; 

5152 + 5253 + S7S8 + 5859 

5154 + S3S6 + 5457 + S659 

5155 + S3S5 + S5S9 + S5S7 

5153 + S7S9 - 2S4S6, 
S1S7 + S3S9 — 2S2S8; 

515254 + SeSgSg + S2S3S6 + S4S7S8 

515255 + 558359 + 5253S5 + 55S7S8 
54S5S7 + 53S5S6 + 55S659 + 51S4S5 
5l5253 + 575359 — 2S4S5S6, 
S1S4S7 + 53S659 — 2S25558, 

515256 + 545359 + 5253S4 + 56S7S8 
S2S4S7 + S3S6S8 + S2S6S9 + S1S4S8 

515257 + S3S8S9 + S2S3S9 + S1S7S8 

545759 + S1S3SG + SeSrSg + S1S3S4 

5153S5 + 55S7S9 — 545058 — 5254S6, 
51S5S7 + 53S5S9 — 5254^8 — S25658; 

51525457 + S3S65859 + 5253S6S9 + S1S4S7S8 — SiS25558 — S255S859 — 5253S5S8 
54575359 + 5iS25356 + 5657S3S9 + S1S25354 - 54S5S657 - 53S4S5S6 - 54S5S659 
515253S5 + S5S7S8S9 - S4S5S658 - 5254S5S6, 

S1S4S5S7 + S3S5S6S9 - S2S4S5S8 - S2S5S6S8, 

51525458 + S2S6S8S9 + S2S3S6S8 + S2S4S7S8 - S1S2S5S7 - S3S5S8S9 - S2S3S5S9 
S4S657S3 + S253S4S6 + 54S65859 + 5iS254S6 - 54S5S7S9 - S1S3S5S6 - S5S65759 
525355S3S9 + S1S25557S8 - S2535653S9 - S1S254S7S3, 
S4555657S9 + S1S3S45556 - 54S6575859 - 5iS2535456, 

525355S658 + 5254S5S7S8 + 5iS25455S8 + 5255S65859 

— S3S5S65859 — 51S2545557 — S1S4S5S7S8 — 5253S5S659, 

5455565359 + 51S254S5S6 + 54S5S657S8 + 5253S4S5S6 

- 51S253S5S6 - 54S5S7S859 - S1S2535455 - S5S6575859. 

In Section 4, it is sufficient to consider only first nine zero-potentials, i.e. we put a. 
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S1S2S8 

545657 



555657 - 51S5S6 

S2S5S9 — S3S5S8 
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S25557S8, 
S1S45556, 



S1S5S7S8, 
S1S3S4S5; 
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for i ^ 10. 

In the case of n.n. and n.n.n. hopping anisotropy (considered in Section 5), the 
zero-potentials invariant with respect to 7r-rotations are of the form: 

Here the zero-potentials, in the case considered in Section 5, are presented: 



,(10 

(11 
T 

(12 
T 

(13: 

T 

l(i4; 

l(15 
l(16 
(17 



(is; 

T 

(19; 

T 

(2o; 



k 



T 
(21 
T 
, (22 

, (23 

rC'p 
l(24 

l(25 

1.(26: 

, (27; 

fVrp 

, (28; 

A/fp 

l(29 

ri/rj-i 



S4S5 - 

S2S5 - 

S1S5 - 
2S4S6, 
2s2S8; 
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=S4 -|- ~ 2S5; 

=S2S3 + S7S8 
-S\S^ + SgSg 
= S3S6 + S4S7 
= S2S4 SeSs 

=52*6 + 
=SlS3 S7S9 
--S]_S-j S3S9 
=SlS2S5 + S5S8S9 
= S2S5-S6 + S4S5S8 
=SiS4S5 + S5S6S9 
= S2S3-S5 + S5S7S8 
-SzS^Sq S4S5S7 
=S2S4S5 + 

=S2S3S4 + S6S7S8 
=S2S4S7 + S3S6S8 
=SlS4S8 + S2S6S9 
=SlS2S3 + 5758^9 
=SiS4S7 S3S6S9 
=SlS2S7 -|- S3S8S9 
=SiS3S6 + S4S7S9 

=SiS2S8 + S2S%Soi 
=SlS3S5 + S5S7S9 
=SlS5"S7 + S3S5S9 



•55 

S5S8, 
S5S7, 



S2S3S6 - 


S4S7S8, 


S2S3S6 - 


S4S7S8, 


■S2S3S6 - 


S4S7S8, 


S1S2S4 — 


•S6S8S9, 


S1S2S4 — 


S6S8S9, 


S1S2S4 — 


S6S8S9, 


S1S2S6 - 


S4S8S9, 


S3S4S5 - 


S5S6S7, 


S2S5S7 - 
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S2S5S9 - 


S1S5S8, 


2S4S5S6, 




2S2S5S8, 




S2S3S8 ~ 


•52*758) 


S1S4S6 — 


S4S6S9, 


S1S3S4 - 


S6S7'59) 


S1S7S8 - 


=52=53^9, 


S2S4S6 ~ 


S4S6S8, 


"S2S4S8 ~ 


"S2'S6'S8; 
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fVrp 


— S1S2S4S5 + S^SqSsSq 


— S2S3S5S6 


- 8485878s 


^(32) 

fVrj-1 


^81828487 + 838G8s8g 


— S2S3S5S8 


— 8285878s 


^{33) 


=SiS4SYSg + S2S3S6S9 


— S1S2S5S8 


— 525558-59 


1.(34) 


= SlS2S3S6 + ■S4S7S8S9 


— S1S4S5S6 


— S4S5S6S9 


^(35) 


=81828384 + SgSTSsSg 


— 83848585 


— 54555657 


7.(36) 


=SlS2S3S5 + S5S7S8S9 


— S2S4S5S6 


— 54S556S8 


2.(37) 


=SlS4S5S7 + S3S5S6S9 


— S2S4S5S8 


— S2S5S6S8 


fbrp 


=SiS2S4S8 + 828q8s8q 


— S1S5S7S8 


— S2S355S9 


jL(39) 


=82838q8s + S2S4S7S8 


— S1S2S5S7 


— S3S5S8S9 


1.(40) 

A/p> 


=S2S3S4S6 + 84858788 


— 81838485 


— 55555759 


l(41) 


= SiS2S4Sq + SiSQSgSg 


— S1S3S5SQ 


— S4S5S7S9 



7 (^4Z_) 

Krp —5152-535456 + 5450575359 — 5153545555 — 5455555759, 
7 (43) _ , 

Krp — 5i52555758 i" 5253555359 — 5i52545758 — 5253565859, 

7 (44) _ . 

—5152545557 + 5355555359 — 5253555553 — 5254555753, 

k^^^ =5152545553 + 5255555359 - 5154555753 - 5253555559, 
; (46) , 

/C^ =5152535555 + 5455575359 - 5i52545555 - 5455555359, 
; (47) 

=5152535455 + 5555575359 - 5253545555 - 5455555753. 

In Section 5, it is sufficient to consider only first nine zero-potentials, i.e. we put ctj = 
for I ^ 10. 
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Appendix B 



1 (102) 2 (101) 3 (100) 4 (99) 5 (98) 6 (97) 7 (96) 8 (95) 9 ^94^; 10 (93) 



xnn 



XX 



11 (92) 



12 (91) 13 fPO; 14 15 (88) 16 f^T) 17 (86) 18 f55; 19 20 (83) 



21 f52; 22 (81) 23 24 (79) 25 (^75; 26 (^77; 27 (76) 28 ^75; 29 (74) 30 (^73; 



S 2J H3 H3 H3 H=: 



31 (72) 32 (^77; 33 ^70; 34 (69) 35 36 (67) 37 (^dd; 38 (^dj; 39 (64) 40 ^^ij 



41 (62) 42 (^di; 43 ^^O; 44 (59) 45 f55; 46 (57) 47 f5di) 48 ^55; 49 (54) 50 fJi; 

• I * 



• I » 

51 (52) 



lire 29 ! All the T-plaquette configurations, over which the potentials II84I and 1851 are minimized, up to the 
symmetries of H (in the isotropic case of n.n. hopping only). The configurations that can be obtained from the displayed 
ones by the hole-particle transformation are not shown, only the numbers assigned to them are given in brackets. 
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63(106) 64(105) 65(104) 66(103) 67(102) 68(101) 69(100) 70(99) 
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73(96) 



74(95) 75(94) 76(93) 



81(88) 



82(87) 



77(92) 
» I » 



78(91) 79(90) 80(89) 



83(86) 84(85) 



U.r6 30 i All the T-plaquette configurations in the case of anisotropic interactions of n.n. hopping, up to rotations 
by TT and reflections in lattice lines parallel to the axes. The configurations that can be obtained from the displayed ones 
by the hole-particle transformation are not shown, only the numbers assigned to them are given in the brackets. 
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Fig U.r6 31 i All the admissible T-plaquette configurations, up to rotations by tt, that are used in constructing the 
fourth-order phase diagram in a neighborhood of point a of Fig. I2UI i.e. the elements of 




1 2 3 4 5 6 



Fie UrG 32; All the admissible T-plaquette configurations that are used in constructing the fourth-order phase diagram, 
up to rotations by tt, in a neighborhood of point b of Fig. 1211 i.e. the elements of (the vertical n.n. pairs that are 
occupied by one ion are forbidden). 
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Fie U.r6 33 i All the admissible T-plaquette configurations, up to rotations by tt and reflections in lattice lines parallel 
to the axes, that are used in constructing the fourth-order phase diagram in a neighborhood of point c~ of Fig. 1211 i.e. the 
elements of (the horizontal n.n. pairs that are occupied by two ions are forbidden). 
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Appendix C 



In this section, at each point p = {u!,e) we give the sets {ai{p)} of the zero-potential 
coefficients a, that we used to construct the four phase diagrams, presented in the pre- 
vious section. Specifically, for each of the four phase diagrams we define a partition, 
called a-partition, of the (u;,£:)-plane (see Figs. OIHSZI displayed below and Tables IHHTTl in 
Appendix E) into sets that differ, in general, from the domains §£>. To each set of such a 
partition we assign a set {(oj, bi, q)} of triplets of rational numbers, such that the coeffi- 
cients ai{p) (with p in the considered set) are affine functions of the form aiUJ+biS+Ci. We 
have not been able to find one set {(a^, hi, Cj)}, that turns the fourth order potentials into 
m-potentials in the whole (cj, £:)-plane. We have not also succeeded in assigning to each 
domain S^i an exactly one set {(a^, 6j, q)}. The same remarks apply to the (cj, 5)-plane. 
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Figure 34: The sets A, B, C, whose boundaries are marked with dashed hnes, used in Table|Hl 
to define the a-partition of the fourth order phase diagram in the case of fermions and for 
5 = O.The dashed-hne segments are determined by their intersection points: a = (—1/2,-6), 
b = (-13/2,-10), c = (19/2,-6), d = (31/2,-10). The coordinates of the remaining points 
can be read from Fig. 111! 
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Figure 35: The sets A, B, C, whose boundaries are marked with dashed hnes, used in Table|Hl 
to define the a-partition of the fourth order phase diagram in the case of hardcore bosons and 
for 5 = 0. The dashed-line segments are determined by their intersection points: a = (1/2,4), 
b = (2, 1), c = (-5/2, -2), d = (9/2, 4), e = (3, 1), f = (15/2, -2). 
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Figure 36: The sets A,B,. . . , whose boundaries are marked with dashed and continuous hues 
(the continuous hnes are the boundaries of the phase diagram), used in Table Uni to define the 
a-partition of the fourth order phase diagram in the case of fermions and for e = 0. 




Figure 37: The sets A,B,. . . , whose boundaries are marked with dashed and continuous hnes 
(the continuous hnes are the boundaries of the phase diagram), used in Table ITTI to define the 
a-partition of the fourth order phase diagram in the case of hardcore bosons and for e = 0. 
The dashed-line segments are determined by their intersection points: a = (—5/2,-1), b = 
(3/2, -3), c = (7/2, -13/3), d = (9/2, -5), e = (6, -4), f = (29/4, -3). 
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Appendix D 



Here we provide equations of line boundaries between the phase domains, and coordinates 
of the crossing points of the hne boundaries that appear in the phase diagrams presented 
in Section 4. The symbols like §i|S2 stand for the line boundary between the phases Si 
and §2, etc. 

Table 1: Domain boundaries of the phase diagrams for fermions, shown in Fig. El 



Lines 


< /32 ^ 7 


7 ^ /32 ^ 10 


^seg 1 ^cb 


= A + 1 




e = 2a; - 2/^2 + 3 




Sidd\S>cb 


e = -2uj + 2p2 + 21 






e = w + i 








w = 2/32 - 1 






- T 




S^2 


e = 2/32 - 6 




CD 1 CD 

^diI^d2 


e = w-| 


CD 1 CD 
^D2l^dl 


e = -a; + 2/32 + ^ 


CD 1 CD 




= /32 + 1 



Points 


< /32 ^ 7 


7 ^ /?2 ^ 10 


a 


(/?2 + i,12) 


(/32 + |,/32 + 5) 


b 


(2/32 - 1,2/32 - 2) 


(/32 + |,/32 + l) 


c 


(f, 2/32 -2) 




d 


(2/32 - 1,2/32 - 6) 




e 


(f, 2/32 -6) 





78 



Table 2: Domain boundaries of the phase diagrams for hardcore bosons, shown in Fig. 1161 



Lines 


< /32 ^ 2 


2 < /32 ^ 5 




w = /32 + f 


^seg 1 Spcfe 


e = 2u; - 2/^2 + 4 




Spcfe Scb 


e = -2uj + 2/32 + 14 






e = w + | 




e = -(J + 2^2 + f 




CJ = 2/^2 + ^ 






w = 1 






e = 2/52 + 1 




QV 1 CD 


e = w + i 


CD 1 CD 
^D2l^dl 


e = -a; + 2^2 + ^ 


CD 1 CD 
^Dll^dl 




w = /32 + 1 



Points 


< /32 ^ 2 


2 < /32 ^ 5 


a 


(/32 + |,9) 


(/32 + |,/32 + 7) 


b 


(2/32 + ^,2/32 + 5) 


(/32 + |,/?2+3) 


c 


(1,2/32 + 5) 




d 


(2/32 + i, 2/32 + 1) 




e 


(1,2/32 + 1) 
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Table 3: Domain boundaries of the phase diagrams shown in Fig. 



Lines 


< /32 ^ 1 


1< ft ^ 2 


2 < /32 5? 3 


3 < ft ^ 6 


S+|S- 


5 = 


Sdd 1 S- 


X _ 1 /32 3 
" ~ 2^2 4 




Sdd 1 §1 


i _ 1^ I /32 11 
6 6 12 




Sdd 1 §2 






Sdd 1 §3 


X — 02 21 
" ~ 2 2 4 




S_|§1 


5 = + ft - 5 


Sl|§2 


5 = -t.; + ft + 2 


S2|§3 


5 = -a; + ft + 6 


§3 1 Sc6 


5 = ~uj + f32 + ^ 






" 4 8 


Sd2 1 ^dd 




i _ " 1 /32 9 
" ~ 4 2 8 




Sd2lS- 




r _ 3u 1 02 3 
" ~ 4 2 8 


SS2IS1 




r _ 3/32 I 1 
" ~ 4, 2 8 


Sdil^a 




I _ 1 27 
" ~ 4 2 "I" 8 


SJIll Sdd 




<^ = 2ft - f 




S^ilSdd 




, , — 23 




Sdl|Sc6 




a; = 2/32 + 


QV 1 QV 




"48 


^vsl^dd 




" 2^2 4 




S^3lS- 




■5 = -2 + 4 


s3;3iSd2 




5= id _ Mi . 17 

" 4 2^8 


QV 1 QV 

^d4l^d2 




X _ 3w I /32 I 13 
" ~ 4 2 8 


Sd4|Sdd 




"22 








<5 = -ft + 5 


QV 1 QV 

^d4l*d3 




i 3w I /32 1 37 

" ~ 4 2 8 


SrfslSdd 




X _ ^ I ;32 29 
" ~ 2 2 4 




SS3IS3 




5 = f - 2ft - 1 


QV 1 

^dsl^dl 




I _ 3cj 1 /32 _i_ 57 
" ~ 4 "I" 2 8 


Gi; 1 QV 




X w 1 7 
'^=-2 +4 


QV 1 QV 

^u2l^d4 




.5 = -ft + 2 


QV 1 Ci; 

^u2l^d3 




5 = f - 2ft + i 


QV 1 ClJ 




a; = 2ft + ^ 
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Points 


< /32 ^ 1 


1< /32 ^ 2 


2 < /32 ^ 3 


3 < /32 ^ 6 


a 


(/32-i,o) 


(-1^0) 


b 




(2/32 - f , ~Y + 1) 


(2/?2-i,f-i) 


(l>o) 


c 




('23 A I n 

V 2 ' 2 2/ 


(23 132 3\ 
V 2 ' 2 2/ 


(i-l) 


d 


(/^2 + f,-i) 


(2/32+f,-/?2 + l) 


V2' 2J 


(2/32+^,-/32 + 3) 


e 


(/32 + f,0) 


V2' 2 2/ 


(4/32-f,-f + 1) 


(2/32-1,-/32 + 2) 


f 




(4/32-|,-f-l) 


(2/32 - 1,-/32 + 1) 


(2/32 + 1,-/32 + 2) 


g 




(2/32-1,-/32) 


(2/32 + f,-/32 + |) 


(2/32 - 1,-/32 + 1) 


h 






ni /32 5\ 

V 2 ' 2 2/ 


(2/32 + f,-/32 + |) 


i 




(/^2+f,-i) 


(19 p2 5\ 
V 2 ' 2 2/ 


(2/32 + f,-/32 + l) 


j 




(2/32 + f , -P2 + 1) 


(2/32 + |,-/32 + i) 


k 




(2/32 + 1,-/32 + 1) 


(2/32 + ^,-/32 + 1) 


1 




(2/32 + ^,-/32 + i) 


(2/32- I, -/?2) 


m 




(2/32-i,-/?2) 





81 



Table 4: Domain boundaries 



Lines 


U < P2 ^ o 


Q 1 Q 

(D_ 1 0^ 


X n 


e 1 Qv 


A a; 9 


QV 1 QV 


- 4 + 8 


QV 1 QV 


X a; 1 
= - 2 - 4 


C 1 Cf 


X LU 7 


<?V 1 QV 


r u 3/32 23 
" — A 2 8 


QV 1 QV 


X StJ I p2 \ b 

^ = -—+2+S 


^-\^d2 


X 3lj 1 /32 11 


QV 1 Q 


r 3/32 39 
" — A 2 8 


QV 1 QV 

^v2\^d3 


r o /3 15 
= 2 - 2/^2 - ^ 


CD 1 Q 


A O 25 


QV 1 Qi; 


r 3w 1 /32 1 25 


QV 1 Q 


r w /32 5 
" ~ 4 2 8 


Cl) 1 QV 

^v2\^dl 


LJ — Zp2 + "2" 




w = 2/32 + f 


S_ |§i 


5 = -a; + /32 - i 


§1 §3 


(5 = + /32 + 2 


§3 Scfe 


6=-LJ + p2 + f 



phase diagrams shown in Fig. UHl 



Points 


U < /52 ^ o 


a 


(-2'0) 


D 


( 1 n\ 
l-2'OJ 


C 


f_5 _1) 
\ 2' 2/ 


d 


(2/32 + ^,-/32-1) 


e 


(2/32 + 1,-/32-2) 


f 


(2/32 + 1,-/32-1) 


g 


(2/32 + f,-/32-|) 


h 


(2/32 + f ,-/32-3) 


i 


(2/32 + ^,-/32-7/2) 


j 


(2/32 + 1,-/32-4) 



Table 5: Domain boundaries of the phase diagrams shown in Fig. |221and Fig. I2HI 



Lines 


fermions 


bosons 


^seg 1 Scfe 


w = 1 


a; = f 


s;;iis:^i 


w = 1 


^ = 1 




e = uj + l 


e = w + | 




e = -co + f 




e = uj-l 


e = u + \ 




e = -uJ + ^■ 



Points 


fermions 


bosons 


a 


(1,5) 


(1,7) 


b 


(1,1) 


(1,3) 



Table 6: Domain boundaries of the phase diagrams shown in Fig. |2l]and Fig. 12^1 



Lines 


fermions 


bosons 








§_|S+ 


5 -- 


= 




Points 


fermions 


bosons 




" 4 8 


" 4 8 




a 


(-io) 


(-1,0) 


s:;iis^3 


"48 


b = - + - 

" 4^8 




b 


(1,0) 


(-1,0) 


s:;3is- 


Si 1 1 


5 — -- -- 

" 2 4 




c 


^2' 2) 


2' 2) 




" 2^4 


" 2 4 
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Table 7: Domain boundaries of the phase diagrams shown in Fig. 1211 and Fig. |2ZI 



Lines 


fermions 


bosons 


s:;3is- 


6 — -- + - 


6 — -- - ^ 

" 2 4 


QV 1 QV 

^v3\^v2 


" 2^4 


6 = - ^ 

" 2 4 


QV 1 QV 

^v3\^d2 


" 4^8 


e w 23 
" ~ 4 8 




j: 3u) 3 
" ~ 4 8 


r 3a; 11 
" ~ 4 8 


QV 1 QV 


j: I 13 

" ~ 4 8 





QV 1 QV 

^d2\^d3 




" 4^8 




" 4^8 


r w 39 
" ~ 4 8 


QV 1 QV 

"dA\"v2 


6 = 2 




Sd4|S2 


6 = 1 




QV 1 Gil 


r 3uj , 37 
" ~ 4 8 




Qt) 1 QV 

^d3\^v2 


6 — - + - 


A — _ 15 
2 4 


Sd3l§3 


"24 


e w 25 
" ~ 2 4 


Qt) 1 QV 

^d3\^dl 


j: 3(i; I 57 
" ~ 4 8 


s: 3a; j 25 
" ~ 4 8 




A - 1 27 
" 4 8 


6 — -- -- 

" 4 8 


Gt) 1 QV 

^dl\^v2 








§_|8i 




§l|S2 


6 = -uj + 2 




§1 83 




6 = -uj + 2 


§2 §3 


S = -uj + 6 




§3 Scfe 


6 = ~co + f 





Points 


fermions 


bosons 


a 






b 




(1,-2) 


c 




fs _5^ 

V2' 2/ 


d 


\ 2' 2/ 




e 






f 




fii -1) 

V 2 ' 2/ 


g 


f3 1) 

V2' 2/ 




h 


(11 1) 
I 2 ' 2/ 




i 


(-io) 
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Appendix E 



Table 8: The case of isotropic n.n. hopping only. The set of zero-potential coefficients 
{oij} in the case of fermions and for 5 = 0. In the first column the sets of the a-partition 
are specified. For more comments see the text in Appendix C. 











0(4. 




A 








oj 1 3 
96 64 





1 

48 










5a; £ 3 
192 64 128 


a) 1 £ 1 9 
32 32 64 


a) 1 £ 1 19 
64 64 384 


S,inB 








192 ^ (IT ^ T28 


- + ' 1 
32 ' :-!2 Gi 


- + ^ 11 
()1 ^ ()1 384 


S^snB 








5tj 1 e 1 21 
192 ^ 64 ^ 128 


uj e 15 
32 32 64 


OJ £ 53 
64 64 384 










5a; _l_ e _l_ 33 
192 64 128 


o) _|_ £ I 27 
32 32 64 


t<j _l_ £ _l_ 73 
64 64 384 


§„2nc 








5w £ _|_ 9 
192 64 ~'~ 128 


UJ £ I 3 

32 32 ' 64 


u> £ _|_ 1 
64 64 ~'~ 384 


Sdinc 








U! 1 6 13 

192 64 128 


O) 1 £ 19 
32 32 64 


U! 1 £ 65 

64 64 384 



Tabic 9: The case of isotropic n.n. hopping only. The set of zero-potential coefficients 
{ai} in the case of hardcore bosons and for 6 = 0. In the first column the sets of the 
a-partitioii arc Ki)ecifie(l. Lor iiicn'c coiiiiiiciitK sec the text in Ap],)cii(lix C. 







«2 


as 




as 


A 








U! 1 5 

96 192 





1 

48 


SpcftHB 








UJ £ 1 23 
24 64 192 


o) 1 7£ 37 
16 ' 192 192 


o) 1 5£ 91 
32 ' 256 768 


S^iHB 








u) £ 1 37 
192 64 384 


7uj _|_ 7£ 55 
192 192 384 


5o> _|_ 5£ 143 
256 ~'~ 256 1536 


S^2nB 








24 ~ 32 ~ 96 


OJ £ 3 

16 16 32 


OJ £ 13 

32 32 192 










_|_ £ _|_ 17 

24 ~'~ 64 ~'~ 192 


o) 1 7£ _|_ 23 
16 ~'~ 192 ~'~ 192 


o) 1 5£ 1 29 
32 256 ~'~ 768 


§.2nc 








OJ £ _|_ 13 

24 32 96 


OJ 

16 16 32 


OJ £ 1 17 

32 32 192 


Sdinc 








w , £ 47 
192 64 384 


7oj I 7£ 125 
192 192 384 


5a) 1 5£ 293 
256 ~'~ 256 1536 
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Tabic 10: The case of isotropic n.n. hopping only. The set of zero-potential coefficients 
{qi} in the case of fermions and for e = 0. In the first column the sets of the a-partition 



are specified. For more comments see the text in Appendix C. 







02 


as 


04 


as 


S_nD 


3200 




^9(T ^ 0400 ^ 04 





1 

4S 


S_nE 


w 5 I 3 
32 16 64 


(5 
8 


uj 1 3 
96 64 





1 

48 


S_nF 


Lj 5,5 
48 16 96 


(5 
8 


1 

24 





1 

48 


SinG 


Uuj 9S , 121 
960 160 ^ 1920 


(5 

8 


7ui . 6 . 403 
5760 320 11520 


UJ S 11 

ISO 80 360 


UJ 1 1 

144 288 


SinH 


3ui 13(5 _|_ 11 
160 160 160 


40 10 ^ 80 


13w (5 I 17 
960 320 320 


LJ 5 1 
40 40 80 


1 

48 


§201 


9S , 21 
128 256 


6(5 I 9 
64 128 


w , (5 .19 
96 128 256 


35 17 
128 256 


u 13 
612 3072 


§20 J 


UJ 7(5 1 3 
128 128 64 


oi 9(5 1 3 
64 64 32 


UJ 1 1 
96 16 


u _l_ (5 5 
128 128 64 


UJ 1 (5 7 
128 128 192 


§3 


7uj 17(5 1 49 
384 192 ' 256 


(5 
8 


29cj 7(5 1 107 
2304 1152 ^ 1536 


144 ^ 144 96 


w 5 1 5 
144 144 ' 96 


Sc6 


359(5 
5728 


(5 
8 


173u) 1 3(5 1 495 
17184 ' 11456 ' 11456 





1 

48 




i5 
8 


(5 
8 


UJ , S , '6 
96 96 64 


(5 

48 


(5 1 
48 48 


SddnB 


.5 
8 


(5 
8 


96 ^ 24 


(5 

48 


(5 1 
48 48 


SddHC 


.5 
8 


6 
8 


Ul 6:3 

96 96 64 


5 

48 


6 1 
48 48 
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Table 11: The case of isotropic n.n. hopping only. The set of zero-potential coefficients 
{cKj} in the case of hardcore bosons and for e = 0. In the first column the sets of the 
a-partition are specified. The cases, where the set §d\t is a proper subset of §>td{p) are 



marked by the asterisk. For more comments see the text in Appendix C. 





ai 


a2 


"3 


0:4 


as 


S_ni 


"16 


§ 


UJ 1 

~96 ~ 192 





~48 


S-nJ 


uj <5 1 9 
132 ~ 16 ^ 64 


s 

~ 8 


7uj ( 3(5 _|_ 91 
192 32 384 


3uj 3(5 27 
~ 64 ~ 32 ^ 128 


1 

^48 


S_nK 


3uj 5 7 
~T28 ~ 16 ~ 256 


3lj 5 23 
64 ~ 8 + 128 


(jj 1 (5 1 19 
48 + T(i + 9fi 


UJ 3(5 21 
32 64 128 


15uj 7(5 331 
512 256 3072 


S_nL 


uj S 5 
16 ~ 16 ~ 32 


^ 5,7 
"■§■ ~ 8 16 


lluj 3(5 203 
384 + 32 ' 7(i8 


2luj . 5(5 I 13 
128" 32 + 256 


41uj 19,5 1027 
~'5T2 ~ 12s ~ 3072 


S_nM 


<5 1 1 
■~T6 + 16 


5(5 1 
~32 ~ 32 ~ 64 


23uj 5(5 1 25 
~ 768 ~ 256 + 1536 


7lj 1 7(5 1 7 
128 128 + 256 


UJ S 19 
"128 ~ 128 ~ 768 


SiPlS 


9oj lid 189 
160 40 320 


UJ 1 (5 1 21 
20 20 40 


5 65 
192 16 384 


3uj <5 43 
80 10 160 


UJ 35 133 
80 40 480 


SinT 


oj 7S , 1 
40 80 ^ 20 


UJ S 1 1 
40 10 ^ 80 


UJ (5 , 11 
60 160 ^ 480 


UJ S 1 
40 40 80 


1 

48 


SanP 


u} S 5 
192 12 384 


6 
8 


UJ 5.5 83 
288 144 576 


w a 67 
192 24 384 


UJ (5 17 
288 72 576 


§3nQ 


2Bui S 475 
192 6 384 


UJ 5 I 35 
8 16 32 


3uj 7(5 269 

128 192 768 


UJ (5 37 

16 32 64 


5uj (5 33 
96 24 64 




6 

16 


UJ 5(5 I 13 
72 36 144 


"144 288 288 


UJ I (5 13 
72 72 144 


1 

48 




5 

16 


S 
8 


UJ 1 5 
96 192 





1 

48 




36 
16 


6 
8 


UJ 1 37 
192 384 


UJ (5 9 
32 16 64 


UJ (5 35 
64 8 384 




5 

32 


3S 
16 


UJ (5,5 
32 96 64 


1 

8 


1 

12 


S„2nB * 


4 11 
32 16 64 


UJ (5 1 11 
16 8 "'' 32 


24 ^ 96 


1 
8 


1 

12 




S 1 
192 24 384 


UJ (5 1 
32 4 64 


7uj (5 , 29 
192 48 ~'~ 384 


u; (5 25 
96 48 192 


1 

12 


SdinC 


S 

16 


(5 
8 


47 

192 384 


3u) 49 
64 128 


49 

32 192 


SdinD 


6 

16 


(5 
8 


UJ 47 
192 384 


UJ 19 
64 128 


1 

48 


Sd2nN 


cj S 31 
64 8 128 


5uj S , 67 
64 16 ~'~ 128 


24 ^ 96 


S 1 
16 4 


UJ S 133 
128 32 768 


Sd2no 


3u; .5 _|_ 9 
32 8 64 


5uj (5 3 
64 16 128 


UJ I 13 
24 96 


3uj (5 11 
64 16 128 


3uj .5 49 
128 32 768 


Sd3nE * 


w 5 11 
32 16 64 


ui 3(5 1 1 
20 20 4 


UJ S 5 
60 20 96 


UJ (5 23 

32 16 64 


3uj 3(5 43 
160 80 192 


Sd3nF * 


S 

16 


3u; 3(5 1 5 
160 20 64 


31t^ (5 1 13 
960 20 384 


(5 3 
16 16 


UJ 3(5 53 
320 80 384 


Sd3nG * 


3u} 1 ^ 1 3 
64 16 128 


OJ 5,3 
16 4 "f" 32 


OJ 1 5 1 5 
240 "T" 20 48 


3uj , 6 21 
64 16 128 


9uj , 3(5 41 
320 80 384 


Sd3nH * 


UJ 1 5 1 25 
64 16 128 


OJ 5 5 
32 4 64 


19uj 1 (5 1 73 
960 20 384 


UJ 1 (5 1 1 
64 "T" 16 128 


UJ 1 3(5 1 
80 80 48 



87 



Below, in a series of tables, we provide the coefficients {oi}, i = 1,...,9, of the 
zero-potentials for the phase diagrams presented in Section 4. The coefficients that are 
missing in a table are equal to zero. As in Appendix D, the symbol S1IS2 denotes the 
boundary between the phases Si and §2, etc. 

Table 12: Zero-potentials coefficients for the phase diagram of fermions, shown in Fig. 1161 
for (32 = 0. 



{ai} 






§1)2 Sdl 


a\ 





012 










0^4 


u; 1 3 
96 ~^ 64 


96 ~'~ 64 


w 1 
96 64 




uj 1 3 
96 64 


a; _|_ 7 
96 64 


w 1 
96 64 


"6 





1 

8 




1 

48 


1 

12 


as 


1 

48 


1 

12 



Table 13: Zero-potentials coefficients for the phase diagram of fermions, shown in Fig. ll6[ 
for /?2 = 7/2. 





QIC Q 1 Qf C 1 Qv Qv 1 QV QV 1 Cli 
dseglOcfe aseg|Oj,i OcblO^i ^v2\^vl ^v2\^dl 


a\ 





012 





as 





04 


LU 1 3 

96 ~'~ 64 


"5 


1 23 
96 ~^ 192 


"6 





Of 


1 

48 


as 


1 

48 
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Table 14: Zero-potentials coefficients for the phase diagram of fermions, shown in Fig. ll6( 
for (32 = 7. 



{Oii} 


QIC Q 1 QV Q 1 QV 

^seg\"cb "seg\"vl "cbl^dl 


QV 1 QV QV 1 QV 

^i;2l^i;l ^v2\^dl 


ai 





Q!2 





"3 





^4 


UJ 1 

96 64 


uj 1 3 
96 64 


as 


w _l_ 49 
96 ~'~ 192 


1 37 
96 ~'~ 192 




1 
8 





aj 


1 

12 


1 

48 


as 


1 

12 


1 

48 



Table 15: Zero-potentials coefficients for the phase diagram of fermions, shown in Fig. ll6[ 
for /?2 = 10. 



{ai} 


QIC? Q 1 QV Q 1 QV 

^seg\"cb "seg\"vl ^cbl^dl 


QV 1 QV QV 1 QV 
^v2\^vl ^v2\^dl 












as 





a4 


1 

6 


a; _j_ 3 
96 ~'~ 64 


"5 


1 

6 


UJ , 49 
96 ~'~ 192 


"6 


1 
8 





a^ 


1 

12 


1 

48 


as 


1 

12 


1 

48 



89 



Table 16: Zero-potentials coefficients for the phase diagram of hardcore bosons, shown 
in Fig. Uni for (32 = 0- 







§1)2 §ul 


§ t;2 §(il 







a2 










^4 


1 5 
96 ~'~ 192 


w _l_ 17 
96 ~'~ 192 


w 7 
96 192 


"5 


li) 1 5 
96 ~'~ 192 


w _l_ 17 
96 ~'~ 192 


uj 7 
96 192 







1 

8 


ay 


1 

48 


1 

12 


"8 


1 

48 


1 

12 



Table 17: Zero-potentials coefficients for the phase diagram of hardcore bosons, shown 
in Fig. Cni for (32 = L 







Q 1 QV C 1 Qtl QV 1 <?V QV 1 QV 

^seg\v^l Ocb\0^i ^v2\^vl ^v2\^dl 


ai 





a2 





as 





a4 


1 

24 


uj , 5 
96 ~'~ 192 


"5 


1 

24 


uj 1 3 
96 ' 64 


"6 





ay 


1 

48 


as 


1 

48 



90 



Table 18: Zero-potentials coefficients for the phase diagram of hardcore bosons, shown 
in Fig. eg for (32 = 2. 



{Oii} 






QV 1 QV QV 1 QV 

^i;2l^i;l ^v2\^dl 


ai 





Q!2 





"3 





^4 


1 

12 


U) 1 

96 192 


iU 1 1 

96 ~'~ 24 


as 


1 

12 


oj _l_ 19 
96 ~'~ 192 




1 

8 


1 

32 


1 

32 


aj 


1 

12 


1 

48 


1 
192 


as 


1 

12 


1 

48 


5 

96 



Table 19: Zero-potentials coefficients for the phase diagram of hardcore bosons, shown 
in Fig. Uni for j32 = 5. 



{ai} 


QIC? Q 1 QV Q 1 QV 

^seg\"cb "seg\"vl ^cbl^dl 


QV 1 QV QV 1 QV 
^v2\^vl ^v2\^dl 












as 





a4 


11 

96 


a; _|_ 1 
96 ~'~ 24 


"5 


11 
96 


1 19 
96 ~'~ 192 


"6 


1 
8 


1 

32 


a^ 


1 

12 


1 
192 


as 


1 

12 


5 

96 



91 



Table 20: Zero-potentials coefficients for the phase diagram shown in Fig. El for j32 





§+|S_ 


S_ 1 §1 


§1 1§2 


§2 1 §3 


§3 Scfe 


5 = 0, ^ f 







uj 1 3 
16 ~'~ 32 


uj 3 
16 32 


ui 9 
16 32 


uj 21 
16 32 










UJ 1 1 
8 ~'~ 16 


uj 3 
8 16 


uj 3 
8 4 


uj 21 
8 16 





as 





uj 1 1 
8 16 


ll) 3 
8 16 


a; 3 
8 4 


u) 21 
8 16 





0:4 


w 1 3 
96 64 


li) 1 1 
96 16 


96 64 


"5 


96 64 




li) 1 1 
96 16 


96 64 


"6 





1 

16 


1 

32 





aj 


1 

48 


1 

96 


1 

48 


as 


1 

48 


1 

96 


1 

48 



Table 21: Zero-potentials coefficients for the phase diagram shown in Fig. El for h 





S-i- 1 §— 


§_|Sl 


Sl|§2 


§2 §3 


§3 Sc6 


5 = 0,uj^f 


ai 





UI 1 1 
16 ~'~ 32 


UJ 3 
16 32 


UJ 11 
16 32 


UJ 23 
16 32 










UJ 1 

8 16 


UJ 11 
8 32 


UJ 41 
8 48 


UJ 23 
8 16 





as 





UJ 1 

8 16 


ui 13 
8 32 


UJ 43 
8 48 


UJ 23 
8 16 







LU , 3 

96 64 


UJ , 9 
96 ~'~ 128 


UI , 13 
96 ~'~ 192 


a; _j_ 3 
96 64 


as 


oj 1 13 
96 192 


UJ , 29 
96 ' 384 


UI I 5 
96 ~'~ 64 


_j_ 13 

96 ~'~ 192 


ae 





1 

32 





a-j 


1 

48 


1 

96 


1 
192 


1 

48 


as 


1 

48 


1 

96 


5 
192 


1 

48 



92 



Table 22: Zero-potentials coefficients for the phase diagram shown in Fig.^l for (32 = 3/2. 



lines & points 


ai 


a2 


as 


Q;4 


"5 


ae 


a7 


as 


"9 


S+|§_ 











w 1 3 
96 ^ 64 


UJ , 5 
96 64 


1 

64 


7 
192 


1 

48 





S_|Si 


16 


Ul 1 

8 8 


UJ 3 
8 32 


1 15 
96 256 


UJ i b 
96 64 


1 

256 


17 

384 


1 

192 





§1 1§2 


49 

16 256 


UJ 31 
8 64 


UJ 13 
8 32 


UJ , 7 
96 128 


OJ 1 27 
96 256 


7 

256 


23 
384 


23 
384 


19 
256 


§2|§3 


uj 105 
16 256 


UJ 115 
8 128 


UJ 29 
8 32 


UJ , 7 
96 128 


w 1 43 
96 512 


25 
512 


157 
3072 


113 
3072 





§3 §c6 


ui 3 
16 4 


UJ 3 
8 2 


a; 3 
8 2 


UJ , 3 
96 64 


UJ , 5 
96 64 





1 

48 


1 

48 
















I 3 
96 64 


UJ , d 
96 64 


1 

64 


1 

48 


7 
192 


















1 

24 


7 

96 





1 

48 


1 

48 


















7 

96 


1 

24 





1 

48 


1 

48 





b 


1 

32 





1 

16 


1 

24 


7 

96 





1 

48 


1 

48 





c 


1 

32 


1 

16 





7 

96 


1 

24 





1 

48 


1 

48 





e 


3 

32 


1 

16 


1 

8 


1 

32 


1 

16 





1 

48 


1 

48 





f 


7 

64 


1 

8 


3 

16 


1 

96 


11 
192 


1 

64 


5 

96 


5 

96 


3 

64 



Table 23: Zero-potentials coefficients for the phase diagram shown in Fig. El for = 2. 



lines & points 


ai 


a2 


as 


Q!4 


a5 


ae 


Q!7 


as 


ag 


S+|§_ 











UJ 1 3 
96 ' 64 


w 1 17 
96 ~'~ 192 


1 

32 


5 

96 


1 

48 





S_|§i 


UJ 1 

16 32 


UJ 3 
8 16 


UJ 3 
8 16 


UJ 1 3 
96 64 


w 1 17 
96 ' 192 





1 

48 


1 

48 





Sl|§2 


UJ 1 

16 4 


UJ 9 
8 16 


UJ 7 
8 16 


UJ , 3 
96 64 


UJ , 23 
96 ~'~ 192 


1 

32 


1 

12 


1 

12 


3 

32 


S2IS3 


UJ 15 
16 32 


UJ 15 
8 16 


UJ 15 
8 16 


UJ 1 3 
96 64 


UJ , 17 
96 ~'~ 192 


1 

16 


1 

12 


1 

24 





§3 §cb 


UJ 25 
16 32 


UJ 25 
8 16 


UJ 25 
8 16 


UJ 1 3 
96 ' 64 


UJ , 17 
96 ~'~ 192 





1 

48 


1 

48 





(5 = 0, f 











UJ I 3 
96 ' 64 


UJ , 17 
96 192 


1 

32 


1 

48 


5 

96 
















1 

32 


7 

96 





1 

48 


1 

48 





(f.o) 











7 

96 


1 

32 





1 

48 


1 

48 





b(e) 


1 

16 





1 

8 


1 

32 


7 

96 





1 

48 


1 

48 





c 


1 

16 


1 

8 





7 

96 


1 

32 





1 

48 


1 

48 






93 



Table 24: Zero-potentials coefficients for the phase diagram shown in Fig. El for /?2 = 5/2. 



lines & points 


"1 


a2 


"3 


a4 


"5 




"7 


as 


ag 


§+!§_ 











cj 1 3 
96 ^ 64 


a> 1 19 
96 192 


1 

32 


5 

96 


1 

48 





§_|Si 


bJ 1 

16 16 


UJ 1 

8 4 


Ul 1 

8 4 


UJ _|_ 3 
96 64 


UJ 1 19 
96 ~'~ 192 





1 

48 


1 

48 





§1 [§2 


uj 9 
16 32 


UJ 5 
8 8 


UJ 1 

8 2 


UJ _|_ 3 
96 64 


UJ 1 25 
96 192 


1 

32 


1 

12 


1 

12 


3 

32 


§2|§3 


u) 33 
16 64 


w 31 
8 32 


- - 1 

8 


96 64 


UJ I 35 
96 384 


9 
128 


1 

12 


5 
192 





§3 §c6 


u) 13 
16 16 


ui 13 
8 8 


ui 13 
8 8 


UJ 1 3 
96 64 


UJ 1 19 
96 192 





1 

48 


1 

48 
















UJ 1 3 
96 64 


UJ 1 19 
96 192 


1 

32 


1 

48 


5 

96 


















1 

48 


7 

96 





1 

48 


1 

48 


















7 

96 


1 

48 





1 

48 


1 

48 





b 


1 

16 


1 

4 


1 

4 


5 

96 





3 
32 


1 

12 


1 

48 





c 


1 

32 


1 

16 





7 

96 


1 

48 





1 

48 


1 

48 





d 


1 

16 





1 

8 


1 

32 


1 

12 





1 

48 


1 

48 





e 


1 

32 


1 

8 


3 

16 


1 

96 


1 

16 


1 

16 


5 

96 


5 

96 





f 


1 

8 


1 
8 


1 

8 


7 
192 


7 

96 


1 

64 


1 

96 


1 

48 





g 


3 

32 


3 

16 


1 

16 


1 

12 


3 

64 


1 

64 


1 

48 


1 

48 





h 


1 

16 


1 
8 


1 

4 


1 

96 


7 

96 


1 

16 


1 

12 


1 

12 


1 

16 


i 


1 

16 


1 
4 


1 

4 


5 

96 





3 
32 


1 

12 


1 

48 






94 



Table 25: Zero-potentials coefficients for the phase diagram shown in Fig. El for = 3 



lines & points 




^2 


^3 


0^4 


as 






as 


s+ 


l§- 











96 64 


96 64 


32 


5 

96 


48 




|§i 


ui 3 
16 ^ 32 


8" ~ 16 


8" ~ 16 


~96 64 


uj 1 7 
~96 64 





1 

~48 


1 

~48 


Si 


§2 


U) 11 
16 32 


LU 7 

8 16 


Ul 11 
8 16 


ll) I 3 
96 ~'~ 64 


LxJ 1 3 

96 ~'~ 64 


1 

16 


1 

12 


1 

12 


§2 


§3 


ui 9 
16 16 


1 

8 


ui 17 
8 16 


ll) I 3 
96 ~'~ 64 


96 32 


5 

64 


1 

12 


1 

96 


§3 




uj 27 
16 32 


uj 27 
8 16 


uj 27 
8 16 


LJ 1 3 
96 64 


uj 1 7 
96 64 





1 

48 


1 

48 


6 = 


= 0, ^ f 











1^ 1 3 
96 64 


a; 1 7 
96 64 


1 

32 


1 

48 


5 

96 




0) (b) 











1 

96 


7 

96 





1 

48 


1 

48 


(f >0) (c) 











7 

96 


1 

96 





1 

48 


1 

48 


d 


1 

16 





1 

8 


1 

32 


3 

32 





1 

48 


1 

48 


e(h) 





1 

4 


1 

4 


1 

96 


5 

96 


1 

8 


1 

12 


1 

12 


i 





1 

4 


1 
4 


5 

96 


1 

96 


1 

8 


1 

12 


1 

12 


f 


5 

32 


7 

32 


3 
16 


35 
768 


7 

96 


3 
256 


5 
384 


1 

48 


g 


1 

8 


1 
4 


1 

8 


3 

32 


3 

64 


1 

64 


1 

48 


1 

48 



95 



Table 26: Zero-potentials coefficients for the phase diagram shown in Fig. El for (32 = & 



lines & points 



02 



as 



0:4 



as 



§+|S- 



"96 ~'~ 64 



96 



11 
64 



32 



_5_ 

"96 



§_|Si 



16 



_9_ 

32 



11 
16 



11 



16 



96 



+ 



64 



96 



11 
64 



48 



§1 1 §5 



16 



19 
32 



11 
16 



15 



16 



96 



64 



96 64 



_1_ 

' 12 



§2 I §2 



16 



25 
32 



21 
16 



23 



16 



96 



+ 



64 



96 64 



_3_ 

32 



12 



§3|§ 



cb 



16 



33 
32 



33 
16 



33 



16 



96 



+ 



64 



96 



11 
64 



"48 



6 = 0,00^ f 
b 

c 
d 

e 



96 



16 

3_ 

16 

T_ 

32 

3_ 

16 

n 

32 

_5_ 

16 



1 

8 

3 
8 

_9_ 

16 

5 
8 

_9_ 

16 

1 

2 



+ 

_!_ 

96 

_ 13 
96 

_!_ 

32 

_13 
96 

11 
192 

_ 11 
96 

1 

192 

__5_ 

32 



64 



96 



11 
64 



13 

96 

_J_ 

96 

_5_ 

32 

5 
192 

_5_ 

96 

_J_ 

48 

2_ 

96 

__3_ 

64 



32 
32 
32 



J_ 

"64 

J_ 

"64 

_3_ 

"32 

J_ 

"64 

_l_ 

"64 



"48 

_5_ 

96 

J_ 

"48 

_1_ 

"48 

J_ 

"48 

J_ 

" 12 

J_ 

" 12 

J_ 

"48 

J_ 

"48 



96 



Table 27: Zero-potentials coefficients for the phase diagram shown in Fig. for = 



liTiPc: AT" 'r\rnri1"c; 


\JL \ 


Lx.2 


"3 












S , IS 


n 


n 

\j 


n 


LU 1 5 

96 ~'~ 192 


LU 1 5 

96 ~'~ 192 


n 


1 

48 


1 

48 


S — 1 § 1 


LO 1 3 

16 ~'~ 32 


LU 1 1 

8 16 


LU 1 1 

8 16 


LU 1 5 

96 ' 192 


LxJ 1 5 

96 ~'~ 192 


n 


1 

48 


1 

48 




LO 9 

16 32 


UJ 1 

8 16 


UJ 1 

8 16 


uj 7 
96 192 


LU 7 

96 192 


1 

8 


1 

12 


1 

12 


'-'3 \ '-'CD 


oj 13 
16 32 


oj 13 
8 16 


oj 13 
8 16 


LU I 5 

96 ~'~ 192 


LU 1 5 

96 ~'~ 192 


n 


1 

48 


1 

48 


A - iu> ^ 


n 


n 


n 


LU , 5 

96 ~'~ 192 


uj 1 5 
96 ~'~ 192 


n 


1 

48 


1 

48 


b 


n 


n 


n 


3 

32 


3 
32 


1 

8 


1 

12 


1 

12 




n 


n 


n 


3 
32 


3 

32 


1 

8 


1 

12 


1 

12 




3 

32 


1 

16 


1 

16 


1 

12 


1 

12 


1 

32 


1 

96 


1 

96 


d 


1 

16 


1 
8 


1 
8 


3 
32 


3 

32 


1 

8 


1 

12 


1 

12 


e 


1 

16 


3 
8 


3 
8 


1 

96 


1 

96 


1 

8 


1 

12 


1 

12 


f 


3 

32 


9 
16 


9 
16 


7 

96 


7 

96 


3 

32 


1 

12 


1 

12 


g 


5 

32 


5 

16 


5 

16 


1 

12 


1 

12 


1 

32 


1 

48 


1 

48 



Table 28: Zero-potentials coefficients for the phase diagram shown in Fig. ^[ for (32 = 3. 



lines & points 




a2 


"3 


04 


"5 


"6 


ai 


"8 


S+|S_ 











a; 1 5 
96 192 


OJ 1 17 
96 ~'~ 192 


1 

32 


5 

96 


1 

48 


§_|Si 


LU 3 

16 32 


LU 5 

8 16 


LU 5 

8 16 


LU 1 5 

96 ~'~ 192 


w 1 17 
96 ~'~ 192 





1 

48 


1 

48 


§1 §3 


LU 15 

16 32 


LU 7 

8 16 


LU 7 

8 16 


ou 7 
96 192 


LU 1 5 

96 ~'~ 192 


1 
8 


1 

12 


1 

12 


§3 §cfe 


LU 19 

16 32 


LU 19 

8 16 


LU 19 

8 16 


LU 1 5 

96 ~'~ 192 


a; 1 17 
96 ~'~ 192 





1 

48 


1 

48 













oj 1 5 
96 ' 192 


OJ _|_ 17 
96 ~'~ 192 


1 

32 


1 

48 


5 

96 


b 











1 

32 


3 

32 


1 

32 


5 

96 


1 

48 













3 

32 


1 

32 


1 

32 


1 

48 


5 

96 


c 


3 

32 


1 

16 


1 

16 


13 
192 


25 
192 


1 

64 


1 

96 


1 

48 


d 


1 

4 


1 

2 


1 

2 


3 

32 


3 

64 


1 

64 


1 

48 


1 

48 


e 


1 

4 


3 
4 


3 
4 


7 

96 


1 

96 


1 

8 


1 

12 


1 

12 


f 


13 
32 


23 
32 


11 
16 


5 
384 


1 

96 


1 

128 


1 

48 


1 

48 


g 


11 
32 


11 
16 


11 
16 


7 

48 


13 
192 


1 

64 


1 

48 


1 

48 



97 



Table 29: Zero-potentials coefficients for the phase diagram shown in Fig. 1241 



lines 


ai 


§+|S_ 







LU 1 1 

16 ~'~ 32 


QV 1 QV 
^v3\^v2 


UJ 7 
16 32 


6 = 0, 00^ I 






Table 30: Zero-potentials coefficients for the phase diagram shown in Fig. 1251 



lines 




§+|S_ 





§-|s;;3 


t<j 1 9 
16 ~'~ 32 


QV 1 QV 
^v3\^v2 


UJ 1 1 

16 ~'~ 32 


5 = 0, ^ -i 






Table 31: Zero-potentials coefficients for the phase diagram shown in Fig. 1261 



lines & points 




02 


"3 


a4 


"5 




^7 


as 


§-|S^3 








UJ 1 1 

8 ~'~ 16 


UJ 1 1 
24 ~'~ 16 


oj 1 3 
96 64 


1 

32 


1 

48 


1 

48 


§_|Si 





1 

8 


UJ 1 1 
4 "T 4 


7<^ 3 
96 64 


96 ~'~ 64 





1 

48 


1 

48 


§1 1§2 


7uj 17 
96 64 


7uj 1 
48 32 


5uj 3 
48 32 





UJ 1 

96 64 


1 

8 


1 

12 


1 

12 


§2 §3 


7uj 29 
96 64 


7uj 21 
48 32 


5uj 19 
48 32 





UJ 1 1 

96 64 


3 
32 


1 

12 


1 

48 


§3 §c6 








UJ 21 
4 8 


7uj _|_ 45 
96 ~'~ 64 


UJ 1 3 
96 ~'~ 64 





1 

48 


1 

48 


§dl Sc6 








5 
4 


6 5 
16 96 


5 

96 





1 

48 


1 

48 


QV 1 QV 
^v2\^dl 








5 
4 


6 1 
16 96 


5 
192 


1 

64 


1 

48 


1 

48 


K3\K2 


1 

8 





UJ 3 
8 16 


UJ 1 3 
24 ~'~ 32 


UJ 1 3 
96 ~'~ 64 


1 
8 


1 

12 


1 

12 


C 


1 

8 





1 
4 


7 

96 


1 

48 


3 
32 


1 

12 


1 

48 


e 








3 
8 


1 

32 


3 

64 


1 

64 


1 

48 


1 

48 



98 



Table 32: Zero-potentials coefficients for the phase diagram shown in Fig. 1271 



lines & points 




^2 


"3 


04 


"5 


"6 




"8 







1 

4 


ui , 9 
8 16 


7 

24 48 


_j_ 5 
96 ~'~ 192 


1 

32 


1 

48 


1 

48 


§_|Si 





1 

8 


LU 1 1 

4 "T 4 


7w 13 
96 192 


LU 1 5 

96 ~'~ 192 





1 

48 


1 

48 


§1 §3 


1 

4 





UJ 1 

4 8 


7uj 1 
96 192 


UJ 7 
96 192 


1 

8 


1 

12 


1 

12 


§3 Scb 








uj 13 
4 8 


7a; 1 83 
96 ~'~ 192 


ui , 5 
96 ~'~ 192 





1 

48 


1 

48 










5 
4 


S 7 
16 96 


7 

96 





1 

48 


1 

48 










<5 
4 


S 1 
16 32 


3 

64 


1 

64 


1 

48 


1 

48 


§:;3is:;2 


1 

8 





8 16 


UJ 5 
24 96 


tt) _l_ 5 
96 ~'~ 192 


1 

8 


1 

12 


1 

12 


d 








5 
8 


23 
96 


13 
192 


1 

64 


1 

48 


1 

48 



99 



Here we present the zero-potential coefficients in the generating points of phase 
diagram shown in Fig. EHl 



Table 33: Zero-potentials coefficients for the phase diagram shown in Fig. 1281 



Points {di, ^2) 


(0,0) 


(2,0) 


(0,1) 


(2,1) 


(0,4) 


(2,4) 


ai 


83 


259 


79 


253 


67 


217 


1536 


3072 


1536 


3072 


1536 


3072 




25 


245 


29 


251 


41 


287 


~ 1536 


3072 


1536 


3072 


1536 


^ 3072 






001 






01 




1536 


3072 


1536 


3072 


1536 


3072 


0:4 


25 


245 


29 


251 


41 


287 


1536 


3072 


1536 


3072 


1536 


3072 




53 
3072 


137 
8192 


19 
1024 


139 
8192 


3 
128 


91 
8192 


"6 


53 


335 


19 


365 


3 


91 


3072 


24576 


1024 


24576 


128 


8192 


0:7 


53 
3072 


401 
8192 


19 
1024 


403 
8192 


3 
128 


1117 
24576 


"8 


53 


571 


19 


589 


3 


139 


3072 


24576 


1024 


24576 


128 


6144 


"9 


19 


101 


19 


13 


19 


1451 


3072 


2048 


3072 


256 


3072 


24576 



100 



